AD-A235  877 


US  Army  Corps 
of  Engineers 

Hydrologic  Engineering  Center 


HEC-4 

Monthly  Streamflow  Simulation 


User’s  Manual 


February  1971 


i'.iivY' 


Hi,***' 


H I 

i  4: 


-"V1 


TE 
:Q91 


O  i  • J 


jio 


--.Vi 
SJj  fe> 

J® 

I 

iyv5. 


Si 


¥ 


Approved  for  Public  Release.  Distribution  Unlimited. 


CPD-4 


HEC-4 


Monthly  Streamflow  Simulation 

User’s  Manual 


February  1971 


Hydrologic  Engineering  Center 
US  Army  Corps  of  Engineers 
609  Second  Street 
Davis,  CA  95616-4687 


(916)  756-1104 


CPD-4 


co  ro  -*■  tocDvjocn^WM 


HEC-4 


Paragraph 

1 


4 

5 

6 

7 

8 


MONTHLY  STREAMFLOW  SIMULATION 

HYDROLOGIC  ENGINEERING  CENTER 
COMPUTER  PROGRAM  723-X6-L2340 


CONTENTS 


Paoe 


ORIGIN  OF  PROGRAM  1 

PURPOSE  OF  PROGRAM  1 

DESCRIPTION  OF  EQUIPMENT  1 

METHODS  OF  COMPUTATION  1 

INPUT 
OUTPUT 

OPERATING  INSTRUCTIONS 
DEFINITIONS  OF  TERMS 
PROPCSED  FUTURE  DEVELOPMENT  8 


EXHIBITS 

DETAILED  EXPLANATION  OF  COMPUTER  PROGRAM 

DESCRIPTION  OF  CROUT'S  METHOD 

INPUT  EXAMPLE 

OUTPUT  EXAMPLE 

DEFINITIONS 

SOURCE  PROGRAM 

INPUT  DATA 

SUMMARY  OF  REQUIRED  CARDS 


1/85 


h'-  CO  CO  CO 


hec-4 


MONTHLY  STREAMFLOW  SIMULATION 

HYDROLOGIC  ENGINEERING  CENTER 
COMPUTER  PROGRAM  723-X6-L2340 


1.  ORIGIN  OF  PROGRAM 

This  program  was  prepared  in  The  Hydrologic  Engineering  Center,  Corps 
of  Engineers.  Up-to-date  information  and  copies  of  source  statement  cards 
for  various  types  of  computers  can  be  obtained  from  the  Center  upon  request 
by  Government  and  cooperating  agencies.  Programs  are  furnished  by  the 
Government  and  are  ac"epted  and  used  by  the  recipient  upon  u'ne  express  unuer- 
standing  that  the  United  States  Government  makes  no  warranties,  express  or 
implied,  concerning  the  accuracy,  completeness,  reliability,  usability,  or 
suitability  for  any  particular  purpose  of  the  information  and  data  contained 
in  the  programs  or  furnished  in  connection  therewith,  and  the  United  States 
shall  be  under  no  liability  whatsoever  to  any  person  by  reason  of  any  use 
made  thereof. 

The  programs  belong  to  the  Government.  Therefore,  the  recipient  further 
agrees  not  to  assert  any  proprietary  rights  therein  or  to  represent  the 
programs  to  anyone  as  other  than  a  Government  program. 

2.  PURPOSE  OF  PROGRAM 

Thin  program  will  analyze  monthly  streamflows  at  a  number  of  inter¬ 
related  stations  to  determine  their  statistical  characteristics  and  will 
generate  a  sequence  of  hypothetical  streamflows  of  any  desired  length  having 
those  characteristics .  It  will  reconstitute  missing  streamflows  on  the 
basis  of  concurrent  flows  observed  at  other  locations  and  will  obtain  maxi¬ 
mum  and  minimum  quantities  for  each  month  and  for  specified  durations  in 
the  recorded,  reconstituted  and  generated  flows.  It  will  also  use  the 
generalized  simulation  model  for  generating  monthly  streamflows  at  ungaged 
locations  based  on  regional  studies.  There  are  many  options  of  using  the 
program  for  various  related  purposes,  and  it  can  be  used  for  other  variables 
such  as  rainfall,  evaporation,  and  water  requirements,  alone  or  in  combination. 

3.  DESCRIPTION  OF  EQUIPMENT 

This  program  requires  a  FORTRAN  IV  compiler,  a  random  number  generator 
(function  RNGEN  included,  see  exhibit  2),  and  a  fairly  large  memory  (64K  on 
the  CDC  6600).  Provision  is  made  for  use  of  three  scratch  tapes,  J  (for 
punched  output),  8  and  9. 

4.  METHODS  OF  COMPUTATION 

a.  In  the  statistical  analysis  portion  of  this  program,  the  flows 
for  each  calendar  month  at  each  station  are  first  incremented  by  1  percent 
of  their  calendar-month  average  in  order  to  prevent  infinite  negative 


logarithms.  This  increment  is  later  subtracted.  The  mean,  standard  deviation 
and  skew  coefficients  for  each  station  and  calendar  month  are  then  computed. 
This  involves  the  following  equations: 
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in  which: 

X  =  Logarithm  of  incremented  monthly  flow 
Q  =  Monthly  recorded  streamfiow 

q  =  Small  increment  of  flow  used  to  prevent  infinite  logarithms 
__  for  months  of  zero  flow 

X  =  Mean  logarithm  of  incremented  monthly  flows 
N  =  Total  years  of  record 

S  =  Unbiased  estimate  of  population  standard  deviation 

g  =  Unbiased  estimate  of  population  skew  coefficient 

i  =  Month  number 

m  =  Year  number 

b.  For  each  station  and  month  with  incomplete  record,  a  search  is 
made  for  longer  records  among  the  stations  used,  to  find  that  which  will 
contribute  most  toward  increasing  the  reliability  of  the  statistics  com¬ 
puted  from  the  incomplete  record.  The  mean  and  standard  deviation  are 
then  adjusted.  Equation  5  is  used  to  compute  the  equivalent  record  required 
to  obtain  statistics  equally  reliable  to  these  adjusted  statistics  and  is 
the  basis  for  selecting  the  best  record  to  be  used  in  the  adjustment. 
Equations  6  and  7  are  the  adjustment  equations. 
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The  primes  indicate  the  long-period  values  and  those  without  primes 
are  based  on  the  same  short  period  for  both  stations  1  and  2,  and: 

N  =  Length  of  record 
R  =  Linear  correlation  coefficient 

c .  Each  individual  flow  is  then  converted  to  a  normalized  standard 
variate,  using  the  following  approximation  of  the  Pearson  Type  III  distribution 


*i,«  -  <Xi,m  -  V  /  S1 

V  ■  [«%V!|*1)1/3  -1]  **1  /6 


t  =  Pearson  Type  III  standard  deviate 
K  =  Normal  standard  deviate 

d.  After  transforming  the  flows  for  all  months  and  stations  to  normal, 
the  gross  (simple)  correlation  coefficients  R  between  all  pairs  of  stations 
for  each  current  and  preceding  calendar  month  are  computed  by  use  of  the 
following  formula: 


f  r  N  p  N  ?  N 

l.  ,  ,  =  i  1  -  I  1  -  (  I  x.  X  )2/(  £  X  £  Xt  T  )1 

i.i-i  l  L  ^  i,m  i-l,m  m=1  i,m  m=1  i-l,m  J 

(N-1)/(N-2)J  ?  (1 


in  which: 


x  =  X  -  X 
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e.  If  there  are  insufficient  simultaneous  observations  of  any  pair  of 
variables  to  compute  a  required  correlation  coefficient,  that  value  must  be 
estimated.  Each  missing  value  is  estimated  by  examining  its  relationship 
to  related  pairs  of  values  in  the  current  and  preceding  month  by  use 
of  the  follcving  formula  using  i,  j,  and  k  subscripts  to  indicate  variables 
used  in  the  gross  correlation. 


\3  -  vw  i  y*1  -o 
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Since,  in  order  to  be  consistent  with  the  two  related  correlation  coefficients, 
the  correlation  coefficient  must  lie  between  the  limits  given  by  equation  11, 
the  lowest  upper  limit  and  highest  lower  limit  are  established  for  all  related 
pears,  emd  the  average  of  these  two  limits  is  taken  as  the  estimated  correla¬ 
tion  coefficient. 

f.  Monthly  streamflows  missing  from  the  records  of  the  various  stations 
are  estimated  for  all  stations  for  each  month  in  turn.  Accordingly,  when¬ 
ever  a  missing  flow  is  being  reconstituted,  there  always  exists  a  valid 
value  for  all  stations  already  examined  that  month  and  for  all  remaining 
stations  in  either  the  current  or  preceding  month.  For  these  remaining 
stations,  the  current  value  is  selected  where  available;  otherwise  the 
preceding  value  is  used.  In  order  to  reconstitute  the  missing  value,  a 
regression  equation  in  terms  of  normal  standard  variates  is  computed  by 
selecting  required  coefficients  from  the  complete  correlation  matrix  for 
that  monrh  and  solving  by  the  Crout  method  (See  exhibit  1).  The  missing 
value  is  computed  from  this  regression  equation,  introducing  a  random 
component  equal  to  the  nondetermination  of  the  equation,  as  discussed  In 
the  streamflow  generation  procedure. 

g.  It  has  been  found  that  valid  use  of  the  regression  technique  requires 
that  all  correlation  coefficients  agree  with  the  data  that  will  be  substituted 
into  the  equations  and  that  the  correlation  coefficients  be  mutually  con¬ 
sistent.  Inconsistency  in  the  correlation  coefficients  causes  the  dependent 
variable  to  be  over-defined  and  is  evidenced  by  a  determination  coefficient 
greater  than  1.0.  If  this  occurs  (because  of  incomplete  data),  the  inde¬ 
pendent  variable  contributing  least  to  the  correlation  Is  dropped,  and  a 

new  regression  equation  is  computed.  This  process  is  repeated  as  necessary 
until  consistency  is  reached  (which  must  occur  by  the  time  that  only  one 
independent  variable  remains).  In  order  to  make  the  correlation  matrix 
consistent  with  the  data  matrix,  all  affected  correlation  coefficients 
are  recomputed  after  each  estimate  of  missing  data. 

h.  Normal  standard  deviates  are  then  converted  to  flows  by  use  of  the 
following  equations; 
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imposing  the  constraint: 


0 


(13) 
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i.  When  the  set  of  flows  is  complete,  all  correlation  matrices  should 
be  consistent  except  for  truncation  errors  in  the  computer,  since  the  data 
arrays  are  complete.  Any  consistency  of  matrices  obtained  in  this  manner 

or  of  matrices  read  into  the  computer  will  result  in  determination  coefficients 
greater  than  1.0.  If  this  occurs,  consistency  of  each  correlation  matrix  is 
assured  by  first  testing  all  combinations  of  triads  of  correlation  coef¬ 
ficients  in  the  current  and  preceding  month  for  all  calendar  months  using 
equation  11  and  raising  the  lowest  of  the  three  coefficients  to  obtain 
a  consistent  triad.  The  test  of  consistency  of  each  complete  matrix  is 
marie  by  recomputing  the  multiple  correlation  coefficient.  If  this  value 
is  greater  than  1.0,  further  adjustment  is  required.  Such  further  adjust¬ 
ment  is  obtained  by  introducing  a  coefficient,  successively  smaller  by 
0.  on  the  radical  in  equation  11  and  repeating  all  triad  consistency 
tests  until  all  matrices  are  consistent.  If  consistency  is  not  reached, 
coefficients  in  each  inconsistent  matrix  are  moved  toward  the  average 
value  of  all  coefficients  in  that  matrix  until  consistency  is  reached. 

j.  Generation  of  hypothetical  streamflows  is  accomplished  by  computing 
a  regression  equation,  by  the  Crout  method  (described  in  exhibit  l)  for 
each  station  and  month  and  then  computing  streamflows  for  each  station  in 
turn  for  one  month  at  a  time  using  the  following  equation.  This  process 

is  started  with  average  values  (zero  deviation)  for  all  stations  in  the 
first  month  and  discarding  the  first  2  years  of  generated  flows . 


V  ■  0iV  *  0>V  '  •  •  • +  0j-iki.j-i  4  0/1-1, j 
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in  which : 


K  =  Monthly  flow  logarithm,  expressed  as  a  normal  standard 
deviate 

/3  =  Beta  coefficient  computed  from  correlation  matrix 

i  =  Month  number 

j  =  Station  number 

n  =  Number  of  interrelated  stations 

R  =  Multiple  correlation  coefficient 

Z  =  Random  number  from  normal  standard  population 

k.  Maximum,  minimum  and  average  flows  are  obtained  for  the  entire 
period  of  flows  as  recorded  and  for  specified  periods  of  reconstituted 
and  generated  flows  by  routine  search  technique. 

l.  Provision  is  also  included  in  this  program  for  use  of  the  generalized 
model  requiring  only  h  generalized  coefficients  for  each  station  (in  place 

of  i8)  and  one  generalized  correlation  coefficient  (in  place  of  1?)  for  each 
pair  of  stations,  in  addition  to  identification  of  wet  and  dry  seasons  for 
each  station.  These  are  defined  as  follows: 

(1)  The  average  value  of  mean  logarithms  of  flows  for  the  wet 
season  (3  months).  This  value  plus  0.2  is  applied  to  the  middle  month 
and  tl;e  average  minus  0.1  is  applied  to  the  other  2  months. 

(2)  The  average  value  of  mean  logarithms  of  flows  for  the  dry 
season  (3  months).  This  is  applied  to  all  3  dry  months.  Mean  logarithms 
'"or  months  between  dry/  and  wet  seasons  are  interpolated  linearly. 

(3)  The  average  standard  deviation  for  all  12  months.  This  is 
aprlied  to  each  of  the  1?  months. 

(4)  The  average  serial  correlation  coefficient,  for  all  1°  montv.s . 
This  value  minus  -15  (but  not  less  than  zero)  is  applied  to  each  wet-season 
month,  and  the  value  plus  .13  (but  not  more  than  .98)  is  applied  to  each 
dry-season  month.  The  average  value  is  applied  to  all  intermediate  months. 

(5)  The  average  interstation  correlation  coefficient  for  all  12 
months  is  applied  to  each  month  for  that  pair  of  stations . 

m.  Because  of  limitations  in  computer  memory  size  and  because  of 
increasing  change  of  computational  instability  with  larger  matrices,  the 
number  of  stations  usable  simultaneously  in  this  program  has  been  limited 
to  10.  However,  the  program  can  reconstitute  and  generate  streamflows  for 
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any  number  of  stations  in  groups  of  10  or  less.  It  will  ordinarily  be 
desirable  to  include  one  or  more  stations  from  earlier  groups  in  each 
successive  group  in  order  to  preserve  important  correlations.  In 
addition  to  providing  flow  data  for  all  stations,  it  is  necessary  to 
designate  NPASS  and  to  follow  each  group  of  flow  data  with  a  standard- 
ormat  card  with  NSTX  (number  of  stations  in  next,  pass  tnat  were  also 
■iced  in  preceding  passes)  and  station  identification  numbers  for  those 
stations.  These  numbers  must  be  listed  in  the  same  sequence  as  their 
data  were  arranged  in  earlier  passes.  Data  for  the  new  stations  for  the 
new  pass  should  then  be  read.  None  of  these  flows  can  occur  in  a  year 
later  than  the  latest  year  for  which  flow  data  occurred  in  the  first  pans. 

n.  As  soon  as  flows  are  reconstituted  for  any  pass,  they  are  read 
onto  the  flow  tape.  After  statistics  are  computed  from  transformed 
reconstituted  flows,  they  are  read  onto  the  statistics  tape  (after 
identification  of  stations  in  the  nass  for  future  reference).  Final 
regression  equation  data  for  each  pass  are  read  onto  the  same  tape  at  the 
same  time  (for  use  in  generation  later).  For  each  new  pass,  the  flow 

and  statistics  tapes  are  searched  separately  for  data  for  those  stations 
already  used  that  also  occur  in  the  new  pass.  In  order  to  read  and  write 
intermittently  and  alternatively  on  the  same  tapes,  it  is  necessary  to 
keen  trar'k  of  tape  records  so  as  to  assure  that  any  read  statement  does 
not  read  beyond  the  record  mark  and  so  that  new  write  statements  occur 
at  the  end  of  all  previous  write  statements  that  are  to  be  saved. 

o.  Once  that  statistics  are  put  on  tape,  they  are  retained  throughout 
the  reconstitution  and  generation  processes.  Flows,  however,  are  saved 
only  for  the  set  of  data  in  which  they  were  reconstituted  or  generated, 
until  the  last  pass  for  that  set  is  completed.  In  the  generation  process, 
it  is  necessary  to  save  the  last  flow  generated  for  each  station  in  one 
set  for  >>se  as  the  antecedent  flow  in  starting  generation  in  the  next 
set.  These  are  saved  in  the  QSTAP  array  with  subscript  I3TAP. 

S .  INPUT 

Input  is  summarized  m  exhibits  7  anti  8.  Data  are  entered  consecu¬ 
tively  on  each  card  using  a  simple  variety  of  formats  to  simplify  punching 
and  handling  cards.  Computed  and  generated  flows  cannot  be  1, 000, 000  units 
or  larger,  and  consequently  must  be  expressed  in  units  that  cannot  exceed 
this  magnitude.  Units  should  be  indicated  on  one  of  the  3  header  cards. 
Column  1  of  each  card  is  reserved  for  card  identification.  These  are  ignored 
by  the  computer  except  for  the  A  in  column  1  of  the  first  header  card ,  which 
is  used  to  identify  the  first  data  card.  An  example  of  input  is  given  in 
exhibit  3.  Certain  inadequacies  of  data  will  abort  the  job  and  waste  input 
cards  until  the  next  card  with  A  in  column  1  is  reached.  A  card  with  A 
in  column  1  followed  by  blank  cards  causes  the  computer  to  stop. 
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OUTPUT 


u. 

Printed  output  includes  key  input  information  for  job  identification 
and  all  results  of  computations.  Generated  flows  are  put  on  magnetic 
tape,  and  computed  statistics  are  punched  on  cards  in  the  format  usable 
later  b-  the  program.  An  example  of  printed  output  is  given  in  exhibit  ^  . 

7 .  OPERATING  INSTRUCTIONS 

Standard  FORTRAN  IV  instructions  and  random  number  generator  are 
required.  No  sense  switches  are  used. 

8.  DEFINITIONS  OF  TERMS 

Terms  used  in  the  program  are  defined  in  exhibit  5  . 

<■>.  PROPOSED  FUTURE  DEVELOPMENT 

There  are  cases  where  the  model  used  herein  does  not  reproduce  historical 
droughts  with  reasonable  frequency.  Consequently,  the  model  is  under 
•ontinuous  study  and  development.  It  is  requested  that  any  user  who  finds 
an  inadequacy  or  desirable  addition  or  modification  notify  The  Hydrologic 
En-ineering  Center. 


EXHIBIT  1 


DETAILED  EXPLANATION 
OF 

COMPUTER  PROGRAM 


GENERAL 

Much  of  the  program  is  explained  by  comment  cards  and  definitions 
of  variables.  Supplementary  explanation  follows,  referring  to  sections 
identified  with  the  indicated  letter  in  column  2  of  a  comment  card. 

SECTION  A 

Correlation  coefficients,  R,  and  beta  coefficients,  B,  are  in  double 
nrecision  for  matrix  inversion  computation,  in  order  to  minimize  computa¬ 
tional  instability.  Correlation  coefficient,  RA,  as  originally  computed 
and  stored,  may  be  defined  in  single  precision.  For  computers  with  word 
length  smaller  than  bits,  many  other  variables  in  this  program  should 
be  in  double  precision. 

When  dimensions  are  changed,  the  corresponding  variable  (starting 
with  K)  should  he  changed  accordingly,  as  these  are  used  to  prevent 
exceedence  of  dimensions .  If  an  excessive  subscript  is  used ,  the  job 
will  be  dumped  until  a  card  with  A  in  column  1  is  encountered,  at  which 
tine  a  new  job  is  automatically  started.  If  5  blank  cards  (with  an  A 
in  column  1  of  the  first)  are  encountered,  the  run  will  be  terminated. 

Job  specification  cards  are  read  in  this  section. 

SECTION  B 

NSTAX  is  number  of  columns  in  correlation  matrix.  These  consist  of 
NoTA  columns  for  the  current -month  values  and  a  similar  number  for  ante¬ 
cedent-month  values.  NSTAA  is  initial  column  number  for  antecedent-month 
coefficients.  These  are  computed  from  NSTA,  which  is  read  in  if  statistics 
are  to  be  provided,  rather  than  computed  from  raw  data.  If  raw  data  are 
to  be  used,  NSTA  is  defined  in  the  program  later  and  NSTAA  and  NSTA  must  be 
also.  Data  for  each  new  pass  are  processed  after  transferring  back  to 
statement  42.  In  the  multipass  operation,  NSTX  is  the  number  of  stations 
used  f rom  previous  passes  and  NSTXX  is  the  subscript  of  the  first  new 
station  for  the  current  pass .  Station  ideatification  for  the  NSTX  stations 
must  be  in  the  order  in  which  data  for  those  stations  were  originally  used, 
because  search  of  data  and  statistics  on  tape  is  made  in  this  order.  Flows 
for  these  stations  are  read  from  tape  IQTAP,  and  corresponding  statistics 
iron  tape  ISTAT.  Variables  LQTAP  and  LSTAT  are  used  to  keep  track  of  tape 
position  for  subsequent  writing. 

♦Provided  through  the  cooperation  of  the  Texas  Water  Development  Board. 


EXHIBIT  1 


Months  are  identified  consecutively  by  the  variable  M  starting  with 
the  month  preceeding  the  first  year  of  data.  Some  quantities  to  be 
accumulated  are  initialized.  Station  combination  data  are  stored  for 
the  purpose  of  obtaining  maximums  and  minimums  (section  D)  of  weighted 
flow  values  later.  Tandem  stations  are  identified  for  cases  where  a 
check  on  consistency  of  generated  quantities  is  deemed  appropriate. 

Station  identification  numbers  are  set  to  a  large  number  so  they  will 
not  be  undefined.  The  flow  array  is  filled  with  -1  values  to  indicate 
missing  values.  For  each  station  and  calendar  month,  the  total  flow  and 
number  of  recorded  values  are  computed  for  computing  a  flow  increment 
and  other  statistics  later.  The  minimum  flow  for  each  station  month 
is  also  computed  in  order  to  avoid  negative  logarithms  later. 

SECTION  C 

Station  data  can  be  read  in  random  order.  Stations  are  identified 
by  subscript  in  the  order  in  which  data  for  each  station  sure  first  read. 

The  year  subscript  is  computed.  Negative  subscripts  will  occur  if  data 
are  for  years  earlier  than  the  starting  year  indicated  on  B  card,  and 
data  for  these  are  rejected,  with  diagnostic  printout.  The  stations 
are  counted  and  the  flows  for  each  month  at  each  station  are  counted  for 
the  purpose  of  computing  frequency  statistics  later.  If  the  number  of 
stations  or  years  exceeds  its  dimension  limit,  the  job  is  aborted.  The 
number  of  stations  is  permanently  stored  in  the  NSTNP  array  for  later 
identification  in  multipass  operations.  The  remainder  of  this  section 
is  self  explanatory,  except  to  state  that  permanent  identification 
station  numbers  are  given  for  stations  in  combination,  for  tandem  stations, 
and  for  consistency-test  stations,  and  subscripts  are  identified  for  rapid 
computation  later. 

SECTION  D 

In  this  section,  maximum  and  minimum  recorded  flows  for  each  calendar 
month  ,  the  water  year  and  for  durations  of  1,  6,  and  54  months,  and  average 
flows  are  computed  for  each  station  and  combination.  Durations  do  not 
span  a  break  in  any  record.  Quantities  are  rounded  off  and  printed  in 
fixed -point  format. 

SECTION  E 

The  logarithm  transform  of  flows  is  accomplished  here.  Missing 
values  are  indicated  by  an  impossibly  large  number  (the  -1  used  for 
missing  flows  is  a  reasonable  logarithm  and  therefore  cannot  be  used 
for  missing  logarithms).  Before  the  log  transform,  the  average  flow 
for  each  calendar  month  at  each  station  is  computed  and  one  ( constrained 
to  a  minimum  of  0.1  flow  unit)  is  added  to  each  flow.  If  the  minimum 
observed  flow  for  that  station  month  is  negative,  that  absolute  value 
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is  also  added  before  the  transform.  After  the  logarithm  transform, 
frequency  statistics  for  each  calendar  month  and  station  are  computed. 

An  increment  needed  to  convert  the  logarithms  to  an  approximately  normal 
distribution  is  also  computed  as  an  alternative  future  transform.  Log¬ 
arithms  to  the  base  10  are  used  so  that  statistics  are  comparable  to 
other  commonly  used  statistics .  A  variable  IRCON  is  set  to  1  if  any 
missing  values  are  encountered,  so  that  the  flow  reconstitution  routine 
will  be  called  later.  A  variable  INDC  is  set  to  1  if  the  first  approxi¬ 
mation  of  increments  causes  any  one  of  the  skew  coefficients  to  be 
smaller  than  0.1  or  larger  than  0.1.  In  an  optional  routine  that  follows, 
the  increment  for  each  station  and  calendar  month  is  adjusted  individually 
and  iteratively  (up  to  14  trials)  until  skew  is  within  0.1  of  zero. 

Stations  with  less  than  three  years  of  data  for  any  calendar  month 
are  deleted,  since  skew  and  correlation  computations  require  at  least 
three  items  of  data. 

S ACTION  t 


Correlation  matrices  are  computed  here  for  the  purpose  of  adjusting 
frequency  statistics  for  short -record  stations.  All  correlation  coef¬ 
ficients  are  first  set  to  -4.0  in  order  to  identify  those  not  computed 
later  for  lack  of  sufficient  observed  data.  Then  accumulations  of  the 
various  quantities  required  are  computed  for  all  items  above  the  main 
diagonal  in  the  correlation  matrix  for  each  month,  using  all  data  common 
to  the  two  stations  involved.  If  more  than  two  items  of  data  are  avail¬ 
able,  the  correlation  coefficients  are  computed.  Coefficients  for  the 
main  diagonal  are  set  to  1.0,  and  those  below  the  main  diagonal  are  set 
equal  to  their  symmetrical  element.  Coefficients  between  the  current 
and  preceding  month's  values  are  similarly  computed.  These  items  con¬ 
stitute  an  extension  of  the  matrix  to  the  right,  which  doubles  its  size, 
and  the  new  portion  is  not  necessarily  symmetrical.  Similar  complete 
arrays  of  average  values  and  root-mean-square  values  for  only  those  log¬ 
arithms  common  to  each  pair  of  stations  are  found  for  later  use  in 
adjusting  statistics. 


A  search  is  then  made  to  determine  the  station  that  would  be  most 
useful  in  adjusting  statistics  for  station  months  with  incomplete  record, 
and  the  means  and  standard  deviations  are  adjusted  in  accordance  with  the 
following  equations : 


S’  =  Sx  ,  (S'  -  S2)  R‘iS1/S2 


4  "  X1  +  "  X2)  RS1/S2 


where  primes  indicate  long-period  values,  subscripts  are  1  for  the  short- 
record  station  and  2  for  the  long-record  station  and, 
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X  =  mean  logarithm 

S  =  standard  deviation  or'  the  logarithms 
R  =  correlation  coefficient. 

An  optional  check  of  consistency  of  standard  deviations  between 
adjacent  stations  for  the  same  month  is  next  made.  This  is  to  assure 
that  frequency  curves  do  not  cross  within  three  standard  deviations 
from  the  mean.  If  there  is  a  conflict,  the  standard  deviation  of  that 
station  designated  in  the  input  data  as  the  dependent  variable  is  modi¬ 
fied  accordingly.  All  frequency  statistics  are  then  printed  out. 

SECTION  G 

All  flows  are  next  standardized  by  subtracting  the  mean  and  dividing 
by  the  standard  deviation  for  the  month  and  station.  An  approximate 
Pearson  Type  III  transform  is  then  applied  as  follows: 

K  =  6  [(.s  gt  +  l)1^  -  l]  /g  +  g/6 

where : 

K  =  normal  standard  deviate 

t  =  Pearson  Type  III  standard  deviate 

g  =  skew  coefficient 

New  correlation  matrices  are  then  computed,  based  on  the  normalized 
variates  ana  using  the  same  standard  procedures  previously  employed  for 
correlating  logarithms.  The  sign  of  the  correlation  coefficient  is 
preserved,  since  the  coefficient  will  be  used  to  establish  regression 
equations .  Correlation  coefficients  are  set  to  zero  if  the  variance 
of  either  variable  approximates  zero,  since  the  computation  of  the  coef¬ 
ficient  Is  highly  unstable  and  since  its  use  would  be  of  little  value. 

SECTION  H 

For  jobs  where  correlation  data  are  given,  the  portion  of  the 
correlation  matrix  above  the  main  diagonal  for  all  months  and  the 
entire  correlation  matrix  relating  current  and  preceding  month's 
values  are  read,  with  a  different  card  for  each  pair  of  stations. 

Values  for  all  12  months  are  contained  on  one  card,  and  the  two  stations 
involved  are  identified  on  the  same  card.  An  automatic  check  is  made 
to  assure  that  cards  are  in  the  required  order  of  columns  and  rows 
in  the  correlation  matrix.  When  generalized  statistics  are  used,  only 
one  correlation  coefficient  for  the  entire  year  is  read,  but  card  order 
is  the  same.  Symmetrical  elements  below  the  main  diagonal  are  then 
filled  in  and  values  af  1.0  are  placed  in  the  main  diagonal. 
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Frequency  statistics  are  then  read,  1  cards  per  station,  with  12 
monthly  values  and  station  identifications  on  each  card.  A  check  is 
made  of  the  station  order,  to  assure  proper  subscripting.  When  gen¬ 
eralized  statistics  are  used,  only  one  card  per  station  is  read,  and 
this  contains  the  maximum  and  minimum  mean  logarithms  and  the  average 
standard  deviation  for  the  year.  The  months  of  maximum  and  minimum 
mean  logarithms  are  also  read  and  converted  to  corresponding  subscripts . 
These  subscripts  will  differ  from  the  calendar  month  number  if  the  year 
used  in  the  study  does  not  begin  with  January. 

SECTION  I 

This  section  searches  for  each  calendar  month  the  entire  correlation 
matrix  to  be  the  right  of  the  main  diagonal  for  misiing  correlation  coef¬ 
ficients  due  to  the  nonexistence  of  at  least  three  years  of  simultaneous 
data  for  the  month.  As  soon  as  a  coefficient  between  two  variables  is 
identified  as  missing,  a  search  of  the  correlation  matrix  is  made  to  find 
established  correlation  coefficients  between  each  of  these  variables  (i 
and  ,1 )  and  any  other  variable  (k).  The  range  within  which  correlation 
between  the  two  variables  must  lie  in  order  to  be  mathematically  consis¬ 
tent  with  the  correlation  with  the  third  variable  is  established  by  use 
of  the  following  equati  m: 

Rij  *  \.i  -V(1V)  <i-\/ 

As  each  successive  third  variable  with  established  correlation  coefficients 
is  found,  the  upper  limit  of  is  constrained  to  the  lowest  of  all  upper 
limits  computed,  and  the  lower  limit  is  constrained  to  the  highest  of  all 
such  lower  limits  computed.  When  the  entire  matrix  has  been  searched  the 
correlation  coefficient  is  estimated  as  the  average  of  these  two  constrained 
limits.  If  this  element  is  above  the  main  diagonal,  the  value  is  also 
entered  for  the  element  symmetrically  across  the  main  diagonal.  The 
search  for  further  missing  correlation  coefficients  is  then  continued. 

SECTION  J 

Where  a  correlation  matrix  is  not  to  be  used  for  reconstituting  data 
but  might  be  inconsistent,  a  triad  consistency  test  can  be  made  in  this 
section.  This  Is  done  by  examining  all  groups  of  three  related  correlation 
coefficients,  and  testing  the  lowest  one  to  determine  whether  it  is  above 
minimum  constraint  established  by  the  equation  in  the  preceding  station. 

If  not,  it  is  raised  to  that  minimum.  When  this  is  done,  it  is  possible 
that  the  adjusted  coefficient  had  already  been  used  in  another  triad  test, 
and  consequently  that  previous  test  would  need  to  be  repeated.  In  order 
to  do  this  properly,  the  entire  matrix  is  searched  up  to  12  NSTA  times, 
where  NSTA  is  the  number  of  stations,  until  a  complete  search  reveals  no 
inconsistent  triad  (INDC  =  0). 
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A  coefficient  FAC  of  the  radical  in  the  equation  is  used  in  order 
to  obtain  complete  matrix  consistency  in  difficult  cases,  whenever 
possible  by  this  means .  A  test  for  overall  consistency  is  made  in 
section  K,  and  if  this  fails.  FAC  is  successively  reduced  by  0.2 
until  overall  consistency  is  reached. 

SECTION  K 

The  test  for  overall  consistency  01  the  correlation  matrix  for 
each  month  is  made  by  constructing  foi  each  station  the  correlation 
matrix  that  would  be  used  in  flow  generation  for  that  station  and  com¬ 
puting  the  multiple  determination  coefficient.  If  the  determination 
coefficient  of  the  matrix  for  any  station  and  any  month  exceeds  1.0,  all 
correlation  matrices  must  be  reexamined,  since  some  coefficients  are 
common  to  two  or  more  matrices.  This  is  done  by  reducing  FAC  in  the 
triad  test  (section  J)  by  0.2  and  repeating  all  triad  tests.  If  FAC  is 
reduced  to  zero  and  consistency  is  not  obtained,  an  index  of  NCB  is  set 
to  1  and  an  averaging  routine  is  used  for  each  inconsistent  matrix.  A 
quantity  SUM  is  computed  as  the  average  of  all  correlation  coefficients 
in  that  matrix,  and  each  element  is  modified  by  multiplying  SUM  by  the 
excess  of  determination  coefficient  and  adding  this  product  to  the 
product  of  the  complement  of  this  multiplier  and  the  value  of  the  ele¬ 
ment  in  the  inconsistent  matrix.  The  averaged  or  smoothed  values  are 
re claced  in  the  complete  matrix  for  the  month,  and  this  requires  some 
careful  manipulation  of  subscripts .  A  new  computation  of  determination 
coefficient  is  made  and  the  smoothing  process  is  repeated  up  to  nine 
times  until  consistency  prevails.  If  this  does  not  occur,  the  job  is 
terminated.  When  consistency  is  established  all  complete  matrices 
are  printed  out  and  essential  elements  are  punched  if  desired. 

SECTION  L 

In  reconstituting  missing  data,  a  search  is  made  for  each  month  of 
record  starting  with  the  first  for  stations  that  have  no  record  during 
that  month  (Q=T).  When  one  is  found,  a  search  of  sill  other  stations  is 
made  to  determine  whether  recorded  or  previously  reconstituted  flows  exist 
for  the  current  month  or,  if  not,  for  the  preceding  month.  If  one  is 
found,  it  will  constitute  an  independent  variable  for  estimating  the 
missing  value,  and  its  value  and  pertinent  correlation  coefficients  are 
stored  in  new  arrays  for  computation  purposes.  The  correlation  coeffi¬ 
cients  with  the  dependent  variable  is  temporarily  stored  in  the  NVAR 
(NSTA-^l)  column  to  assure  that  coefficients  relating  independent  vari¬ 
ables  which  have  sufficient  array  space  (they  cannot  exceed  NSTA  in 
number) .  A  variable  ITEMP  counts  the  number  of  independent  variables 
(stations  for  which  recorded  or  reconstituted  data  are  available).  It 
is  incremented  after  its  set  of  correlation  coefficients  are  stored  in 
the  R  array,  and  is  finally  used  to  relocate  the  correlation  coefficients 
involving  the  dependent  variable.  If  no  independent  variables  with  data 
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are  found,  as  can  happen  in  the  first  month  of  record,  a  correlation  is 
made  with  the  preceding  value  for  the  same  station  and  that  preceding 
value  is  arbitrarily  set  at  the  average  for  the  month.  The  regression 
equation  and  deiermination  coefficient  are  then  computed  using  subroutine 
CROUT.  The  variable  having  the  lowest  absolute  value  of  correlation  with 
the  dependent  variable  is  identified,  and  beta  coefficients  axe  searched 
in  order  to  eliminate  all  unreasonable  coefficients.  In  the  usual  case 
where  the  simple  correlation  coefficient  between  any  variable  and  the 
dependent  variable  is  positive,  unreasonable  coefficients  are  assumed  to 
be  those  larger  than  1.5  or  smaller  than  -.5-  In  the  case  where  the 
variable  correlates  negatively  with  the  dependent  variable ,  the  reason¬ 
able  range  is  -1.5  to  0.5.  If  an  unacceptable  coefficient  is  found,  INDC 
is  set  to  1.  If  this  happens  or  if  the  determination  coefficient  does 
not  lie  between  0  and  1.0,  the  variable  with  the  smallest  correlation 
coefficient  is  eliminated,  the  correlation  array  reconstructed  accordingly, 
and  the  regression  equation  recomputed.  This  process  is  repeated  until 
all  required  conditions  exist.  The  missing  value  is  then  computed  by 
use  of  the  regression  equation  and  adding  a  random  component  normally 
distributed  with  zero  mean  and  with  variance  equal  to  the  error  variance 
of  the  regression  equation. 

As  soon  as  the  missing  value  is  estimated  a  search  is  made  for  all 
established  values  in  the  current  and  preceding  month  with  which  it  is 
to  be  correlated,  and  sums  of  logarithms,  squares,  and  cross  products 
are  incremented  in  preparation  for  recomputing  all  affected  correlation 
coefficients.  After  checking  for  sufficient  (three  years)  record  and 
nonzero  variance,  the  correlation  coefficient  is  recomputed.  If  the 
standard  deviation  of  either  variable  is  very  small,  the  correlation 
coefficient  is  set  to  zero.  If  the  coefficient  is  above  the  main  dia¬ 
gonal  of  the  correlation  matrix,  its  value  is  also  assigned  to  sym¬ 
metrical  element.  Since  estimation  of  a  missing  value  affects  correla¬ 
tion  coefficients  between  variables  in  the  current  and  following  month, 
which  coefficients  are  stored  in  a  different  matrix,  this  process  of 
adjusting  the  correlation  coefficient  is  applied  to  those  values  next. 

SECTION  M 

After  all  flows  are  reconstituted,  the  flow  tape  is  read  until  the 
proper  position  for  writing  the  newly  conmuted  flow  data  on  that  tape  is 
reached,  and  headings  are  printed  for  writing  flews  on  the  printer  later. 
Then  the  standard  deviates  are  converted  to  flows  by  reversing  the  Pearson 
type  III  transform,  multiplying  by  the  standard  deviation,  adding  to  the 
mean  and  taking  the  antilogarithm .  The  increment  is  then  subtracted  and 
if  the  resulting  value  is  negative  for  a  variable  with  zero  lower  limit, 
it  is  set  to  zero.  In  the  case  of  reconstituted  flows,  the  Pearson  Type 
III  transform  is  constrained  so  that  the  excess  of  the  standard  deviate 
over  and  above  2.0  is  multiplied  by  a  maximum  of  O.j  (if  the  standard 


7 


EXHIBIT  1 


diviation  exceeds  0.'-;).  This  simply  prevents  obtaining  unreasonably 
extreme  values  due  to  sampling  errors .  It  is  a  moderation  of  the 
extrapolation  rather  than  an  abrupt  truncation. 

The  test  for  tandem  station  consistency  is  next  made,  and  incon¬ 
sistent  flows  are  identified  for  printout  and  changed  to  the  limit  of 
consistency.  The  downstream  flow  is  inane  consistent  with  the  Siam  of 
upstream  flows.  Flows  are  punched  on  cards,  if  desired,  printed  out, 
and  written  on  the  flow  tape  for  use  in  future  passes.  NQ.TAP  is  incre¬ 
mented  and  represents  the  total  number  of  records  on  the  tape. 

SECTION  N 

After  converting  deviates  to  flows,  the  frequency  statistics  are 
recomputed  in  order  to  agree  accurately  with  observed  end  reconstituted 
data.  If  a  consistency  test  is  called  for,  the  variable  ITRN3  is  set 
to  ?  and  computation  is  transferred  to  near  the  end  of  section  F,  where 
the  test  is  made  and  the  transfer  index  causes  a  return  to  this  portion 
of  the  program.  Adjusted  statistics  are  printed,  and  the  consistent 
correlation  matrix  is  nrinteri  (and,  if  desired,  also  punched)  by  transfer 
to  section  K,  using  ITRII.l  as  a  return  indicator  again.  The  statistics 
are  then  punched,  i1'  desired.  Flows  for  the  snecified  station  combina¬ 
tions  are  then  computed . 

SECTION  0 

Maximum  and  minimum  recorded  flows  are  computed  by  transfer  to 
section  S,  using  ITRNS=1  as  a  return  indicator.  The  variable  ITMP 
keeps  a  record  of  the  remaining  years  whose  maximum  and  minimum  flows 
have  not  been  searched  yet. 

Next,  generalized  statistics  are  computed,  if  desired,  (if  IGNRL 
equals  two).  As  indicated,  straight  averages  of  all  12  monthly  correla¬ 
tion  coefficients  in  every  category  are  taken.  Means  are  averaged  for 
the  three  wettest  consecutive  months  and  the  three  driest  consecutive 
months  and  the  seasonal  timing  noted .  Standard  deviations  for  sill  12 
months  are  averageo .  Generalized  statistics  are  then  printed  out. 

Next,  generalized  statistics  read  in  section  H  are  used  to  compute 
required  arrays  of  statistics.  Skew  and  increments  are  set  to  zero. 

The  mean  for  the  middle  month  of  the  wet  season  is  .2  higher  than  the 
wet  season  average  and  means  for  the  other  two  months  are  .1  lower. 

Means  for  the  dig/  seasons  are  uniform,  and  means  for  the  transition 
seasons  are  interpolated  linearly.  Correlation  coefficients  for  the 
dry  season  are  .IS  higher  (constrained  below  .98)  than  the  annual 
average,  and  those  for  the  wet  season  are  .IS  low  (constrained  above 
zero).  All  of  these  operations  are  in  accord  with  the  generalized 
model  developed  in  HEC. 
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SECTION  P 


After  obtainin'*  monthly  statistics  and  correlation  matrices,  re¬ 
gression  equations  for  each  station  and  calendar  month  are  computed. 

Flows  are  generated  in  the  station  order  in  which  data  or  statistics 
are  read  and  are  generated  for  each  month  at  all  stations  before 
proceeding  to  the  next  month .  Flows  at  each  station  are  correlated 
with  flows  of  the  antecedent  month  at  that  station  and  at  all  stations 
for  which  the  current  month's  flows  have  not  yet  been  generated.  For 
other  stations,  flows  for  the  current  month  are  used. 

Regression  equations  are  computed  in  subroutine  GROUT.  If  any 
correlation  matrix  “ormed  is  inconsistent  (which  should  not  occur  at 
this  stage,  except  for  truncation  of  computated  intermediate  variables,), 
a  transfer  to  section  J  is  effected,  and  consistency  operations  performed 
on  all  correlation  matrices.  After  such  a  transfer,  all  regression 
equations  must  be  recomputed,  since  any  correlation  coefficient  might 
have  changed.  After  this,  only  the  beta  and  alienation  coefficients 
need  be  retained,  in  addition  to  the  frequency  statistics.  To  +he  -”.’lti- 
pass  operation,  these  are  all  written  on  tape  I3T3T  at  this  point. 


SECTION 

A  routine  for  projecting  historical  sequences  into  the  future  is 
employed  here.  Values  of  yPREV  (previous  month's  deviate)  for  each 
station  is  determined  as  the  transform  of  the  flow  for  the  month  pre¬ 
ceding  the  first  month  specified  (by  input  data)  to  be  generated.  The 
variable  MA  is  computed  for  the  subscript  or  0  that  conforms  to  the 
'irst  month  of  projected  flows.  If  the  projected  flow  routine  is  not 
to  ve  used,  the  computer  is  next  set  up  to  generate  two  years  of  flows, 
at  the  end  of  which  synthetic  sequences  will  have  a  virtually  random 
start . 

In  the  multipass  operation,  stations  are  identified  and  all  neces¬ 
sary  statistics  are  contained  in  the  order  needed  on  tape  ISTAT.  In  any 
pass  after  the  first,  "lows  generated  in  earlier  passes  for  the  same 
period  (the  same  sequence  oi’  data)  must  be  read  from  tape  I3TAP,  and  this 
tane  must  he  rewound  before  each  pass  in  order  to  permit  a  complete 
search.  In  any  sequence  after  the  first,  the  preceding  flow  for  the 
"irst  month  to  lie  generated  is  the  last  flow  in  the  preceding  sequence 
or  that  station,  and  these  are  saved  in  the  QSTAP  array  for  multipass 
operation.  If  the  multipass  feature  is  not  used,  all  necessary  statistics 
and  flows  for  generating  are  in  memory'. 


SECTION  R 


In  starting  to  generate 
the  year  number  o"  the  "irst 


flows,  a  variable  JXTMP  is  used  to  identify 
year  or  each  sequence  in  the  multipass 
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operation.  Variables  AVG  and  GDV  are  used  to  compute  the  mean  and 
standard  deviation  of  the  deviates  por  each  flow  sequence.  These  are 
later  used  to  adjust  all  deviates  so  that  the  means  and  standard 
deviations  in  ever.-  generated  sequence  will  be  the  same  as  those  of 
the  historical  sequence. 

Variables  JA  and  IJJ  are  set  up  to  correspond  to  the  first  and 
last  year  of  generation  in  each  successive  sequence,  depending  on 
the  type  of  operation.  MA  has  already  been  set  up  as  the  subscript 
of  Q.  corresponding  to  the  first  month  of  flows  to  be  generated  (for 
use  in  nrojecting  historical  flows  recorded  to  the  current  time). 

CvFREV  for  each  station  has  been  identified  as  the  previous  month’s 
flow  for  that  station.  Flows  are  then  generated  for  each  station, 
using  stored  regression  equations  and  a  random  component.  Each  generated 
flow  is  immediately  entered  into  the  QPRGV  array,  because  its  preceding 
flow  will  never  a~ain  be  used  in  that  pass. 

In  fve  multipass  routine,  flows  (as  deviates)  are  written  on  tape 
at  the  end  of  each  nass ,  and  the  last  dew  for  each  station  is  stored 
in  the  :*3TAF  array  ’"or  use  in  the  next  sequence. 

If  more  than  lp  years  (an  arbitrarily  selected  length)  of  flow 
are  being  generated  in  any  sequence,  deviates  are  adjusted  so  that  their 
mean  is  zero  and  variance  1.0.  Their  unadjusted  mean  and  standard 
deviation  are  printed.  Then  they  are  transformed  to  flews,  and,  if 
called  Tor,  consistency  tests  between  stations  are  made.  For  variables 
with  zero  natural  limit,  a  check  for  negative  values  is  then  made.  Flows 
are  then  printed  and.  if  desired,  punched.  Flow  combinations  are  then 
computed . 

SECTION  o 


Before  computing  maximum  and  minimum  values  of  generated  flows, 
a  oositive  value  of  .IX  is  looked  for  to  assure  that  flows  generated 
are  not  to  be  discarded  (the  first  two  years  generated  for  a  random 
start).  Also,  at  .least  NYMXG  years  must  have  been  generated  before 
maximum  and  minimum  values  are  "omputed  (this  applies  only  when  the 
number  of  years  remaining  for  generation  in  the  last  seauence  does 
not  equal  NYMXG) .  Maximum  sums  are  initiated  at  an  extremely  large 
negative  number  and  minimum  sums  as  an  extremely  large  positive  number 
(T).  Then  a  routine  search  of  Flow  sums  for  the  specified  durations 
at  each  station  is  made  for  the  sequence,  and  results  are  printed  out. 
Since  this  routine  is  used  for  reconstituted  flows  as  well  as  for 
generated  flows,  a  transfer  indicator  is  used  to  determine  whether  the 
next  step  is  bacv  to  the  reconstitution  routine  or  the  generation  routine. 
If-  the  latter,  a  check  is  made  for  the  multipass  routine.  If  all  passes 
are  not  completed,  e  transfer  to  section  is  made.  If  all  passes  are 
completed  for  this  sequence  or  if  the  multipass  routine  is  not  being 
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used ,  a  ''heck  is  made  of  renal  via,  .ears  to  r  e  'one rated .  Ij'  greater 
thai  zero,  a  transfer  to  section  -I  is  made  after  adjusting;  years  yet 
to  ’e  generated .  Otherwise  the  'oh  is  ended  and  a  new  job.  if  any,  is 
started  . 


RANDOM  NUMBER  FUNCTION  RNGLN 


This  random  number  function  is  for  a  binary  machine  and  the  constants  must  be 
computed  according  to  the  number  of  bits  in  an  integer  word.  The  numbers 
generated  are  uniformly  distributed  in  the  interval  0  to  1, 

The  function  is  called  from  the  main  program  bv  a  statement  similar  to  the 
following: 


A  =*  RNGEN  (IX) 

Where  A  is  some  floating  point  variable  name  and  TX  is  some  integer  variable 
name.  The  argument  name  IX  need  not  be  the  same  In  the  main  program  and  the 
function.  The  argument  must  be  initial  iced  to  zero  in  the  main  program.  The 
location  of  the  initializing  statement  is  important  and  depends  on  the  results 
desired.  If  it  is  desired  to  have  different  sets  of  random  numbers  for  each 
of  several  different  sets  of  computations  (jobs)  that  are  run  sequentially  on 
the  same  program,  then  the  argument  must  be  iniLial Lzed  at  the  very  beginning 
of  the  program  and  never  reinitialized.  If  it  is  permissible  to  use  the  same 
sequence  of  random  numbers  for  each  job,  the  argument  must  lie  initialized  at 
the  beginning  of  each  job.  The  advantage  of  this  latter  option  occurs  when  one 
of  the  jobs  must  be  re-run  for  some  minor  reason  as  the  same  random  numbers 
will  be  used  and  the  results  will  be  comparable. 

Three  constants  must  be  computed  bv  tite  following  equations: 

Constant  one  (Cl)  ■*  ^ "  +  3 

p 

Constant  two  (C2)  =  2  ’  -1 

Constant  three  (C3)  =  1./2.'1 

Where:  15  ■  number  of  bits  in  an  integer  word 
The  constants  for  some  of  the  common  computers  are  listed  in  the  following  table: 


CTOirUTER  _  _ 

SIZE  OF 
INTEGER 
WORD 

Cl 

CONS".' ANTS 
.....  ^ 

C3 

GE  200  Series 

19 

1027 

524287 

0.190734863E-05 

GE  400  Series 

23 

4099 

8388607 

0.1 19209 290E-06 

IBM  360  Series 

31 

65539 

2147483647 

0.465661287K-09 

IBM  7040  and 

7090  Series 

35 

262147 

34359738367 

0. 291038 3046E-10 

UNIVAC  1108 

II 

II 

U 

II 

CDC  6000  Series 

48 

16777219 

281474976710655 

11 

0. 3552713678)1-14 
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EXHIBIT  2 
Crout's  Method 


On*  of  the  beet  methods  for  solving  systems  of  linear  equations 
on  desk  calculating  machines  was  developed  by  P.  D.  Crout  in  l?4l. 

This  method  is  based  on  the  elimination  method,  with  the  calculation* 
arranged  in  aystamatio  order  so  as  to  facilitate  their  accomplishment 
on  a  desk  calculator.  In  this  method  the  coefficients  and  constant 
terms  of  the  equations  are  written  in  the  form  of  a  "matrix,"  which  is 
a  rectangular  array  of  quantities  arranged  in  rows  and  column*. 


The  method  is  beet  explained  by  an  example.  Suppose  that  in  a 
multiple  correlation  analysis  it  is  required  to  solve  the  following 
system  of  linear  equations  to  obtain  the  unknown  values  of  b  ,  b,,  b. 
and  b5. 


*2 

b2 

+ 

*2X3 

b3 

2*2*4  b4 

+ 

2*2*5 

b5  “  2xxx2 

£*2x3 

b2 

+ 

*3 

b3 

+ 

Ex  x.  b, 

3  4  4 

+ 

2*3*5 

b5  -  2^*3 

2*2*4 

b2 

+ 

2*3*4 

b3 

+ 

4  \ 

+ 

2*4*5 

b5  -  2*1*4 

2*2x5 

b2 

+ 

2*3*5 

b3 

+ 

2*4*5  b4 

+ 

2*c 

5 

b5  “  Ex1x5 

For  simplicity  let  ua  replace  the  coeffioie^ ' s  of  the  b's  by  the  letter* 
p,  q,  r  and  s(  and  the  oonstant  terms  by  t  etter  t,  using  subscripts 
1,  2,  3  and  4  to  denote  the  respective  equt.  .ones. 


P1  b2  +  4,  b,  +  K  +  »,  bs  -  t 


ll  3 


1  “4 


*1  “5 


p2  b2  ♦  q9  b^  4-  r,  bL  +  »9  b*  -  t 


*2  3 


2  “4 


2  5 


p3  b2  +  q3  b3  +  r3  b4  +  s3  b5 


p4  “2 


b«  + 


44  b3  +  r4  \  *  84  b5 


-  t,. 


A  continuous  check  on  the  computations  a*  they  progress  may  be  obtained 
by  adding  to  the  matrix  of  the  above  system  a  column  of  u's,  such  that 
u-p  +  q  +  r+  s  +  t.  The  matrix  and  check  oolumn  are  written  as 
followst 
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*1  *1 

*2  *2 

*5  U? 

*4  % 

The  eleaenta  p^,  q,,,  r^  and  a^  fora  the  "principal  diagonal"  of 

the  aatrix.  Exaaination  of  the  original  equatlona  shows  that  the  oo- 
effioients  are  ayaaetrloal  about  the  principal  diagonal,  i.e.,  q,  -  p_, 
“  Pj»  *2  m  8i  "  P4»  *2  "  <14*  8j  " 

This  is  characteristic  of  the  ays tea  of  equations  to  be  solved  in  any 
multiple  correlation  analysis.  Because  of  this  nyaaatrjr.  the  ooaputa- 
tions  are  considerably  slaplifled.  While  the  Crout  aathod  aay  be  used 
to  solve  any  ay a tea  of  linear  equations,  the  ooaputatlonal  steps  given 
here  are  applicable  only  to  those  with  symanl 1 1  mil  ooeffioianta. 


The  solution  oonaists  of  two  parts,  vis.,  the  computation  of  a 
"derived  aatrix"  and  the  "back  solution."  Let  the  derived  aatrix  be 
denoted  as  follows i 
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The  general  pattern  of  the  above  computations,  which  may  ba  applied 
to  a  system  containing  any  numbar  of  equations,  is  as  follows! 

(1)  The  first  column  of  the  derived  matrix  is  copied  from  the  first 
oolunn  of  the  given  matrix. 

(2)  The  remaining  elements  in  the  first  row  of  the  derived  matrix 
are  computed  by  dividing  the  corresponding  elements  in  the  first  row  of 
the  given  matrix  by  the  first  element  in  that  row. 

(3)  After  completing  the  nr^  row,  the  remaining  elements  in  the  (n+l)— 
oolumn  ere  computed.  Such  an  element  (X)  equals  the  corresponding  element 
of  the  given  matrix  minus  the  product  of  the  element  immediately  to  the  left 
of  (X)  by  the  element  immediately  above  the  principal  diagonal  in  the  same 
column  as  (X),  minus  the  product  of  the  seoond  element  to  the  left  of  (X) 

by  tbs  seoond  element  above  the  principal  diagonal  in  the  same  oolumn  a*  (X), 
eto.  After  eaoh  el  mien t  below  the  principal  diagonal  is  recorded,  and  while 
that  element  is  still  in  the  oaloulator,  it  is  divided  by  the  element  of 
the  principal  diagonal  whioh  is  in  the  same  oolumn.  The  quotient  is  the 
element  whose  looation  is  symmstrioal  to  (X)  with  respeot  to  tbs  principal 
diagonal. 
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(4)  When  the  elements  in  the  (n+l)-^  column  and  their  symmetrical 
counterparts  have  been  reoorded,  the  (n+l)l£  row  will  be  aomplete  except 
for  the  last  two  elements,  which  are  next  computed.  Such  an  element  (X) 
equals  the  corresponding  element  of  the  given  matrix  minus  the  product 
of  the  element  immediately  above  (x)  by  the  element  imediately  to  the 
left  of  the  principal  diagonal  in  the  earns  row  as  (X),  minus  the  product 
of  the  seoond  element  above  (X)  by  the  eeoond  element  to  the  left  of  the 
principal  diagonal  in  the  same  row  ae  (X),  etc.,  all  divided  by  the 
element  of  the  principal  diagonal  in  the  mm  row  ae  (X). 


The  oheok  column  (U)  of  the  derived  matrix  serves  as  a  continuous 
aha ok  on  the  computations  in  that  each  element  in  the  oolumn  equals 
one  plus  the  sum  of  the  elements  in  the  same  row  to  the  right  of  the 
principal  diagonal.  That  is, 

Ui  -  1  +  Q1  -t-  R1  +  3^^  -t-  T1 


1  +  R„ 


r  S2  .  T2 


1  ♦  s3  +  T3 


U,  -  1  +  T. 


This  oheok  should  be  made  after  completing  each  row. 

The  elements  of  the  derived  matrix  to  the  right  of  the  principal 
diagonal  form  a  system  of  equations  which  may  now  be  used  to  oompute 
the  unknown  values  of  b2,  b^,  b^  and  b^  by  successive  substitution. 

This  ie  known  as  the  "back  solution."  The  computatiops  are  ae  follows: 


b 

to 

b 

b 


5 

4 

5 
2 


S2b5  -  H.b4 

Slb5  ’  V*  '  V} 


It  is  very  important  that  the  computations  be  carried  to  a  suffi¬ 
cient  number  of  digits,  both  in  computing  the  coefficients  and  constant 
terms  of  the  original  equations,  and  in  computing  the  elements  of  the 
derived  matrix.  It  is  possible  for  relatively  small  errors  in  the 
ooeffioients  and  constant  terms  of  the  original  equations  to  result  in 
relatively  large  errors  in  the  computed  solutions  of  the  unknowns.  The 
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greatest  source  of  error  in  computing  the  elements  of  the  derived  matrix 
arises  from  the  loss  of  leading  significant  digits  by  subtraction.  This 
must  be  guarded  against  and  can  be  done  bv  carrying  the  computations  to 
more  figures  than  the  data.  As  a  general  rule,  it  is  recommended  that 
the  coefficients  and  constant  terms  of  the  original  equations  be  carried 
to  a  sufficient  number  of  decimals  to  produce  at  least  five  significant 
digits  in  the  smallest  quantity,  and  that  the  elements  of  the  derived 
matrix  be  carried  to  one  more  decimal  than  this,  but  to  not  less  than 
six  significant  digits. 
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EXHIBIT  5 

DEFINITIONS  -  723-X6-L2340 


AC1 

ACS 

ACS 

ADJ 

ADJ1 

ALCFT(I,K) 

ALOG 

AN  LOG 

ANYRS 

AV(I,K) 

AVG(l,K) 

AVGQ(I) 

AVMN(I) 

AVMX(l) 

B(L) 

BETA(l,K,L) 

BLANK 

CROW 

CSTAC(KX,K) 

DABS 

DQ(I,K) 

DTRMC 

E 

FAC 

I 

IA 

IAKAL 

IENDF 

IGNRL 


IMN(I) 

IMNTH 

IMX(l) 

INDC 

IP 

IPASS 

IPCHft 

IPCHS 

IQ(I) 


-  Alienation  coefficient  for  station  1 

-  Alienation  coefficient  for  station  2 

-  Alienation  coefficient  for  station  3 

-  Plus  sign  indicates  value  smaller  than  upstream  sum  by 
tandem  test 

-  Equal  sign  Indicates  value  adjusted  by  tandem  test 

-  Alienation  coefficient  array 

-  Computer  library  function  of  natural  logarithm 

-  Number  of  logarithms 

-  Number  of  years  of  record 

-  Mean  logarithm 

-  Average  of  the  generated  deviates 

-  Average  monthly  flow  for  a  station 

-  Average  logarithm  of  flow  for  minimum  3  consecutive  months 

-  Average  logarithm  of  flew  for  maximum  3  consecutive  months 

-  Beta  coefficient 

-  Beta  coefficient  for  generation  equation 

-  Blank  space 

-  Program  subroutine  to  solve  simultaneous  equations 

-  Coefficient  by  which  flews  are  multiplied  before  adding  in 
a  combination 

-  Computer  library  function  of  absolute  value  of  double  precision 
number 

-  Increment  of  flow 

-  Determination  coefficient 

-  Letter  E  Indicates  estimated  value 

-  Temporary  factor 

-  Index  for  calendar  month 

-  Indicator  in  column  1  of  first  card  for  each  Job 

-  Indicator,  positive  value  calls  for  analysis 

-  End  of  file  indicator 

-  Indicator,  +  2  calls  for  computing  generalized  statistics 
and  +  1  or  +  2  calls  for  using  generalized  statistics  for 
generating  flews 

-  Month  sequence  number  of  last  month  of  3  driest  consecutive 
months 

-  Calendar  month  number  for  first  month  of  water  year 

•  Month  sequence  number  of  last  month  of  3  wettest  consecutive 
months 

-  Transfer  indicator 

-  Month  number  for  preceding  month 

-  Sequence  number  of  pass  (subset  of  stations) 

-  Indicator,  positive  value  calls  for  writing  discharges  on  tape 

-  Indicator,  positive  value  calls  for  punching  statistics 

-  Fixed-point  conversion  of  flow  values 
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IQ5PAP 

IRCON 

ISKZ 

iot(k,l) 

iota(k) 

iotac(kx,k) 

IOTAH 

ISTAP 

istat 

istn(l) 

iott(k.l) 

istx(l) 

IOTY(L) 

ITKMP 

ITMP 

ITMPP 

ITP 

irons 

IX 

IXX 

IYR 

IYRA 

IYRPJ 

J 

JA 

JTMP(L) 

JTP 

JX 

JXTMP 

K 

KM 

KPASS 

KOTA 

kotac(kx,k) 

KOTAP 

KX 

KYR 

L 

LA 

LQTAP 

LOTAT 

LTMP(L) 

LTP 

LTRA 

LX 

LYHPJ 

M 

MA 


-  Tape  number  tor  storing  flows 

-  Indicator,  positive  value  calls  for  flow  reconstitution 

-  Positive  value  calls  for  varying  flow  increment  (DQ)  to 
mate  skew  zero. 

•  Sequence  manner  of  upstream  station  for  tandem  test 

-  Station  number  " 

-  Station  number  in  a  combination 

-  Temporary  station  number 

-  Station  sequence  number  for  all  passes 

-  Tape  number  for  storing  statistics 

-  Station  number  of  downstream  tandem  statical 

-  Station  number  of  upstream  tandem  station 

-  Station  number  of  Independent  station  for  consistences  test 

-  Station  number  of  dependent  statical  for  consistences  test 

-  Temporary  variable 

-  Temporary  variable 

-  Temporary  variable 

-  Temporary  variable  . 

-  Transfer  indicator 

-  Temporary  variation  of  I 

-  Argument  for  random  number  function 

-  Number  of  current  year 

-  First  year  of  data 

-  Year  of  start  of  flow  projection 

-  Index  for  year 

-  Sequence  number  of  projection  year 

-  Matrix  column  number 

-  Matrix  column  number 

-  Temporary  variation  of  J 

-  Temporary  variation  of  J 

-  Index  far  station 

-  Dimension  limit  for  number  of  consecutive  months 

-  Dimension  limit  for  number  of  passes 

-  Dimension  limit  for  total  number  of  stations 

-  Index  number  of  station  in  a  combination 

-  Dimension  limit  far  total  number  of  stations 

-  Temporary  variation  of  K  or  combination  sequence 

-  Dimension  limit  for  number  of  consecutive  years 

-  Index  for  related  station 

-  Temporary  variation  of  L 

-  Number  of  records  up  to  present  position  on  tape  IQTAP 

-  Number  of  records  up  to  present  position  on  tape  ISTAT 

-  Matrix  row  number 

-  Matrix  row  number 

-  Letter  A 

-  Temporary  variation  of  L 

-  Last  year  of  each  projection 
*  Serial  number  of  month 

-  Sequence  number  of  month  of  projected  flow 
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MO(l ) 

MPASS 

MTHPJ 

MXRCS 

N 

NC 

NCA 

NCAB(l,K,L) 

NCB 

NCOMB 

NCSTY 

NINDP 

NJ 

NLCXJ(I,K) 

NMNMX 

NMXMN 

NPASS 

NPROJ 

NQ 

NQTAP 

NS«(L) 

NSTA 

NSTAA 

NOTAC(KX) 

NSTAT 

N9TAX 

N9TNP( I ) 

NSTX 

NSTXX 

NSUM(K) 

NTNDM 

NVAR 

NYMXG 

NYRG 

NYRS 

Q(K,K) 

QM(I) 

QMIN(I,K) 

QPRBV(I) 

QR(M,K) 

QSTAP(I) 

R(K,L) 

RA(l,K,L) 

RAV(K,L) 


-  Calendar  month  number 

-  Temporary  counter  for  number  of  passes 

-  Calendar  month  of  start  of  each  projection 

-  Number  of  years  in  each  period  for  which  maximum  and  minimum 
recorded  and  reconstituted  flows  are  desired 

-  Serial  number  of  period  of  flows 

-  Counter  to  prevent  continuous  looping 

-  Counter  to  prevent  continuous  looping 

-  Number  of  values  and  cross  products  used  to  compute  correlation 
coefficients 

-  Transfer  indicator 

-  Number  of  combinations  of  stations  max.  and  min.  quantities 
are  to  be  computed 

-  Number  of  consistency  tests 

-  Number  of  independent  variables  in  regression  study 

-  Number  of  years  in  computation  sequence 

-  Number  of  logarithms  used  to  compute  frequency  statistics 

-  Number  of  months  following  dry  season  and  preceding  wet  season 

-  Number  of  months  following  wet  season  and  preceding  dry  season 

-  Total  number  of  passes  in  Job 

-  Number  of  projections  of  future  flows  from  present  conditions 

-  Counter  for  number  of  flows 

-  Total  number  of  records  saved  on  tape  IQTAP 

-  Number  of  upstream  stations  in  tandem  test 

-  Number  of  stations  in  analysis 

-  NSTA  +  1 

-  Number  of  stations  in  a  combination 

-  Total  number  of  records  saved  on  tape  ISTAT 

-  NSTA  +  NSTA 

-  Number  of  stations  in  a  particular  pass 

-  Number  of  stations  in  current  pass  that  occurred  in  preceding 
passes 

-  NSTX  +  1 

-  Number  of  stations  upstream  from  a  station  for  tandem  test 

-  Number  of  tandem  tests 

-  Total  number  of  variable  in  regression  study 

-  Number  of  years  of  generated  flows  in  each  period  for  which 
maximum  and  minimum  flows  are  desired 

-  Total  number  of  years  of  generated  flows 

-  Number  of  years  of  recorded  flows 

-  Monthly  flow 

-  Monthly  flow 

-  Minimum  flow 

-  Flow  for  previous  month 

-  Identification  symbol 

-  Temporary  storage  of  QPREV 

-  Correlation  coefficient  in  a  given  matrix 

-  Correlation  coefficient 

-  Average  correlation  coefficient  for  12  calendar  months 
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RMAX 

RMIN 

RNGEN(IXX) 

R1 

R2 

R3 

SD(I,K) 

SDAV(K) 

SDV(l,K) 

SKEW(I,K) 

SMQ(J,K) 

SQA(I,K,L) 

SQB(I,K,L) 

SUM 

suma(i,k,l) 

sumb(i,k,l) 

T 

TEMP 

TMP 

TMPA 

TMPB 

TMPP 

TP 

X(I) 

XI NCR ( I ) 
XPAB(I,K,L) 


-  Maximum  consistent  correlation  coefficient 

-  Minimum  consistent  correlation  coefficient 

-  Program  random  number  function 

-  Correlation  coefficient  being  tested 

-  Correlation  coefficient  being  tested 

-  Correlation  coefficient  being  tested 

-  Standard  devaition  of  logarithms  for  calendar  month 

-  Average  standard  deviation  for  12  consecutive  months 

-  Standard  deviation  of  the  generated  deviates 

-  Skev  coefficient  of  logarithms  for  calendar  month 

-  Maximum  or  minimum  flow  for  month  or  duration 

-  Sum  of  squares  of  first  variable 

-  Sum  of  squares  of  second  variable 

-  Average  correlation  coefficient  of  matrix 

-  Sum  of  first  variable 

-  Sum  of  second  variable 

-  Large  positive  constant 

-  Temporary  variable 

-  Temporary  variable 

-  Temporary  variable 

-  Temporary  variable 

-  Temporary  variable 

-  Temporary  variable 

-  Value  of  independent  variable  in  regression  equation 

-  Iteration  value  for  flow  increment 

-  Sum  of  cross  products  of  first  and  second  variables 
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ooo 


C  7?3»*6»L2340  monthly  «Tr.tAWFt,'»  StrtJLATIUM  MFC,  C  QF  F»  USA  NOT  t97Q 
C*  *****  *  tI»M*KY  FUNCTIONS  A I  OH  r  DAMS  ************  1008 

*P0GW*H  CRCUT.RNSfN  —  8ft  COMMENTS  tH  WNOEN  1003 

INDEXES  IsCAlEnDa«  MONTH  JaYEAR  wajTA  LRNELAtEO  OTA  HAS'-'CCESalTE  »0*TH  1004 

1000 

DIMENSION  1006 

.ortft),w(to,ti),  iootooo 

.Ai.CFTn2>t0)>Av(l?,!0),*v8()2>10),AveflCt0),Av*N(ia)»AVHK(lC}* 

.BETA (la. 10. lo) ,00(1?, lO),lWN(tQ) ,Tm«(10)»IQ<13) * 
(,I8TA{lO),.ITH».{9),LTMP{lO),N0(12),NCAB<18*tO,20)* 

.NL06(l2.in)*ofl201,10),CH(l*),RMth(t2,lo),QPREV(10),qR(120|*lO)* 

.QBf4P(lOO)#»A(l2,lO*2r),t«AV(lO,t«),30(l2»tO)*SOAVtlO)*SO*H2»t9l» 

,, SKEH(  12. 10), SH« (38, 10)  , BOA ( 1 2, |0,  JO) , SOP (12»  12*20) fOUMAt  12*10*20). 
fSHPB  (  12,  10,?0)  ,  A  (1  0)  *  *1  AC*  (12)  ,XPAS(12,  10,20)  * 

.CffTAC  (2,  1 0, 5),  MT>C  (2,  1C),1.1TN(  10)  ,137  7(10,10),  1977(10)* 
,l3T»(IO),P9TAC(2,tO,5),Aa**(10),hsTAC{2,S),HSTNP(S), 

.  N8t'H  { 1  0, 5)  *NcOpB(5),HTNCH(5),lSTtlO*10fS) 

DOUBLE  PRECISION  R,8  1016 

COMMON  0TRmC,NI»0P*8  1017 

DATA  LT"A/JHa/,BUAMK/1H  /,E/lH£/,*OJ/lM*/,AOjl/lMa/  1010 

10  FORMAT (1K1)  1019 

2o  Format  (i«,is*19x6)  t02oooo 

30  FORMAT (lx, t7, 910)  1021 

•0  F0SMAT(1x,a3,9A«,10AA)  1022 

50  F0RMAT(IX«IJ,I4,12F6.<I)  1023 

60  FORMAT(1*,FT.O,9m«,0)  1021 

re  format  a*.T3,I*,12F6.3>  1023 

*0  FORMAT  ( t v , 1 7 • 12FA.3)  102N 

*0  Format  (tx,(7,12FA.t)  io2( 

IOC  FORM*T(tx.I«,X6.l2!«*I10)  1020 

tic  FORMAT  (Al,A3,9A4,10A«)  1029 

120  FORMAT  (1*. 17,3^8.3,216)  1030000 

130  FORMAT  {/23M  GENERALIZED  3tATl3UCS//l3M  ST1  0T2  RAV)  1031 

140  FORMAT (/3«H  5TA  A  VP*  AVMN  SOA V  MA*HO  MtNMO)  1032 

I3TAT43  .  1034 

IqTAPbi  1035 

KPA3.8R5 

KRTApaiOO  1037 

KSTARIO  1038 

KYRRlOO  1039 

KMaKvR* 12*1  1040 

K3T  Ac(l  10*2 

c  waste  CAROS  UNTIL  AN  A  IN  COLUMN  1,  FINJT  TITLE  CARO  1043 

C  A*  CARO  A  ••  1044 

ISO  REAO (5 , 1 l 0)  IA,(3MO(M,1),M»|,20)  1045 

IF  flA.NE.LTRA)  SO  TO  ISO  1046 

WRITE (6, 10)  1049 

R6A0(S,«n) ((3^0(M,K),M«1,20),N*2,3)  1050 

WRITE  (6,40)  (  (3M(J{M,K)  ,M*l  ,20)  ,4*1,3)  IOSI 

C  •*  CARO  S  CARO  C  **  1032000 

HEAr>(5,  J0)r»HA,lMRTW,rANAL,M*Rc3*NTRG,NYMXG*NPAS3,IPCW9.1PCH8,N3TA  1053 
l,NCOMB.NTNOM,»(CSTY,IGNML,NPRnj,tYKPJ,MTt.P),LYRPj  1035 

C  terminate  WITH  5  PLAN*  CARDS,  an  A  IN  COL  l  CF  FIRST  1056 

ITMPaf AMAL*NVRG  1057 

IF(TTMR.GT.O)G0  TO  ISO  1058 

STOP  1059 

1  AO  WRITE  (A, 170)  NYNS,N3TA,NCnwn* (PASS  1060 

170  FORMAT  t/10M  DIMENSION  EXCEEDED  *3«* 4NNYN8, I4,5X,4mN8TA, 13, SX,  10A1OOO 

15WNC0MR, I3,SX,5WTPA8S, II)  1062000 

80  TO  150  1063 

1«C  WRITE (6, 190)  1064 

19c  Format (/fOSM  itra  imntm  ianal  mxpcs  mtrs  ntmx#  npaj*  ipcno  ipcms  I063 
1  N3TA  NcflrtO  NTN0M  NCST»  7GN«L  NP40J  XTRFJ  MTMPJ  LTRPJ  )  1066 

WRITE  (6.20)  ITS  A,  I.-iMT*-,!  ANAL.  M»RC.«,NrRG,NTM»C,NRA33,l*CH0,  IPCMS,  1087' 

1N9TA, NCCwB.NTMOM.NCSTY, 16**1. NFR0J,ITRPJ*NTWPJ,LYMPJ  1068 

IF  (L»«PJ-I7RA,6E,KTN)  «C  TO  160  1069 

CR  *  *  *  *  *  jeT  CONSTANTS  *••**«♦•*•***#*•••••*  io70 
XVXRO  1071 

MSTAARMSTAaI  1072 

N8TAxRM.nTA*N3TA  1073 

TR99999959.  1074 

IVRARlrRA*!  1075 


EXHI8IT  6 


NST  *»0 

NSTXXaf 

IPA.ASBI 

REN|NO  ISTAT 

N8T*T»0 

UST*T»0 

REN I  NO  TOTAF 

NOTAPafl 

LOTAPaO 

OB  195  JalfKpASS 
HGOPP (J)aO 
HTNOM(j)aO 
l9S  CONTINUE 
SO  TO  270 

SAVE  STATIONS  FROM  PREVIOUS  PASSES  IP  NECESSARY 
200  IPASSBtPASSM 
W»ITE(S,10) 

IP  flPASS.GT.KpASS)  80  TO  ISO 

•  A  CARO  J  ** 

REAr)(5,30)NC0«4»MTH0M,NCSTV,KSTX»  ( 1ST A U) , *•! , NSTX) 

WRITE  (h. 21  fl>  IPASS.  (lSTA(X),Ka|,N8TX) 

21 C  FPRMAT(5H0PAS8  • I3/2SH  STA(S)  FRON  PREVIOUS  PASSES  ,1016) 
NST**«»STX*l 
REWIND  IQTAP 
UOTApao 
REKINO  ISTAT 
HPA»S*l 
RE AO  (ISTAT) 

URTAT»| 

ITP8NVRS*12«1 

ITEnp*NSTNP(nPASS) 

ITMPpao 

DO  250  Xnl.NsTX 

220  REAn(loTiPmHP,{0(N,K),X»l,ITP) 

U0TAP»L0TAP*1 

IF  flSTA(KJ.NE.ITPP)  60  TO  220 
230  ITNPP* ITPPpTt 

IF(ITPPP.CT.1TE  V>  OC  T*  2«0 

READ  (ISTAT)  lTHP,(AVO,N),Sn(t,K),3NEW(l,X)»0Q{t,K>,tBETAU,K,O,L 
l«l.t TEMPI, AECVT(I.K),I«)»>>) 
tSTAT»LSTAT*l 

IF  (ITh*.EO.ISTA(K))  60  TO  250 
GO  TO  230 
240  READ ( ISTAT) 

ESTAT«LSTATTt 
MPASS«MPASS*1 
XTEMPaMSTNP (MPASS) 

ITRPPai) 

GO  TO  23o 
250  Continue 

On  260  Kat.NsTX 
RSIIM  (X ,  IP*S9)aO 
00  260  la  1 , 12 
260  MLOGd.XjaNVpS 
2»0  IF(IANAL.GT.O)  NSTAaNSTX 
On  aso  T»l, 12 
M0(T)aiMKTHTI 
IF(P0(I).LT.l3)G0  TO  260 
Hn(i)»'tp(i)-ia 
2«0  CONTINUE 

|F(NCOMP,te.!»  GO  TO  320 

HCCNB(I»»RS)aNCDM8 

IDENTIFY  STATION  COPSlNATION* 

00  300  Kal.NCOnS 

•*  CARO  0  AA 

REA0(5,30)ITP# (tSTAC(X#L) ,1*1, |TP) 

WRITE  (6,290)  K.ITP,  (ISTAC(X,U),l,Al,ITP) 

290  RORNAT  (/5n  COpB, I2,5H  8TA,|Sl8) 

N8TAC (X , |PAS8)aITP 

.  ....  *•  CARO  e  AA 

READ (5, 60)  TENP.(C8TACCR,I.IP*8S)#1.»I»IYR) 


1134 

1135 

1136 
113T 
1138 
1139000 


EXHXIIT  * 


s* 


.l'*  5 

sgaTTEf!'.  7t' .1  (C3T*C(K,L.IPA3fl),L»!,rrP) 

Fn«H»r  t?» laps, 3) 

U4Rnoo 

IF(NtNDM.LF.O)  GO  TO  330 

1145 

fmOM{rP*ss>=>*T*io* 

00  340  L  Xs  1 «  nTNDM 

1  146 

c 

CARD 

F 

#* 

1147 

REA0(9,3C1  IsTN(LX)  « IT ‘Ip,  f I3TT  (L X ,L) .L«l » IT«Pl 

U46 

WRlTE(6.33o)  I.X.lSTN(L<),(jSTTa^,L)#L*l»rT«P) 

1144 

jic 

FORMAT  </13H  TANQFh  GROUP, 13, 6x,UhOQwNSTR£AH  3Ta.i5.6X, 

1130000 

115MUP3T9FA''  STA(Sl,10I5) 

lismoo 

3*0 

NSM*(LXl=ITHP 

1132 

i^O 

IF(IPASS.E7.1 )NYrt3»0 

1153 

on  380  Ksn.-,TxX,KsTA 

1154 

NSUH(K.TPAsS)»0 

ISTA{K)stf>00-K 

1156 

c 

TfitTIATt  -1,  NO  RECORD  FOR  ALL  FLOPS 

1157 

TO  360  Msl.KH 

1138 

7(M,K)*-1 . 

1154 

On  370 

1160 

«Lrr.fi,K)=o 

1161 

00(I,K)»G. 

1162 

SMI* Jt,X)ST 

1163 

37C 

CONTINUE 

1164 

3*C 

Continue 

1165 

tF(NCSTY.LE.C)  Go  TO  «20 

116b 

MUTE  {8. 3001 

1167 

3«C 

Fn«MATC/30v,flHgTATinNS/17W  ChNSISTEnCY  TEST # 5* , 23* I ROEPENDE N T 

OE 

1166 

»nENOENT) 

1169 

00  400  t.sl,NcSTV 

1170 

c 

*  * 

card 

G 

*# 

1171 

«'E»0  (5.30)  Isrx(L),T3TY(L) 

1172 

'i9ITE(b..'Hol  L,  tSTX(L)iISTy(L) 

1173 

«l  0 

Format (1 3X,  13,8 X, I5.8X, 15) 

llTq 

/l?C 

IF(tANAL.lr.o)GO  T(j  1070 

1175 

tc  * 

♦  *  *  *  *  R ?. A 0  AND  PRQCE33  1  STATI()N“YEAR  OF  OAT*  •  «  «  * 

1  *  • 

« 

*  « 

1176 

c 

** 

CARO 

H 

t* 

1177 

><FAn(5,5o)  T3UNfIY»,(aH(I),lai,\2) 

1178 

c 

** 

CARO 

1 

** 

1174 

c 

SLAnK  C tRO  IMITCATES  E*D  CF  FLOh  D*TA 

lific 

rF(r3TAN.LT.l75n  TO  s^o 

lift 

1F(NSTA.LT.1)G0  Tf  45C 

1182 

c 

A(,OlON  SUBSCRIPT  TO  STATION 

1183 

DO  440  K5N.STXX.N5TA 

1164 

IFJIJTAN.E J.I3TA(K))GC  TO  480 

1185 

«TC 

CnNTINljE 

1186 

a?ri 

NRT  ASN3TA* I 

1167 

IF  (NSTA.nT.K3TA)  GO  TO  160 

IS87-2* 

K*N3Ti 

1106 

ISTAfK)3ISTAH 

1189 

c 

ASSIGN  SUBSCRIPT  TO  YEAR 

1190 

J»IYR.IYM 

1191 

:F(NYRS.LT.3.AN0.IPA5T,e0.1)NTBS»J 

1 1  92 

IF(J.GT.n.AN0.J.LE.NYP37  GO  TO  «80 

1193 

WRITE  (b.«7n)IYS 

1  194 

<i 7  0 

format  (/Iom  UNACCEPTABLE.  YEAR, IS) 

1145000 

'in  to  150 

1196 

c 

STORE  FlON3  IN  3TATTPN  AND  MONTH  ARRAY 

1197 

a»c 

MoJ»l,->.l  ) 

1198 

nn  45o  1=1.12 

1  199 

*a,M  ♦  l 

1200 

IF  (ON  ( I )  .1  F, (”!,))  go  rr  »50 

1201 

IF  (OH  (I)  .t.T.)M|N(I,K))  OIN(t,K)*OlCl) 

1202 

Nuor  (i,Kls‘.Liir.(i,K)Tt 

1203 

on(l.KI=n0(r,F)*cH(n 

1204 

G(H,K)«n-(t) 

1203 

noc 

continue 

1206 

-■r  to  »3n 

1207 

«oc 

n9TAA»n.3IA*1 

1208 

IF  (NYRS.GT ,KYM.OR.HSTA»NC0HB,GT.KftTA)  CO  TO  160 

1209 

IF(N3Ta.lE  .  0  )  r.o  TO  18  0 

1210 

NRTNPf XPAS5)sUSTA 

1211 

NSTAX»N3T»*Nr,TA 

1212 

•3* 

EXHIBIT  6 

L 


IFInCOMB.LE.OIGO  to  s«o 

c  identify  an  subscripts  for  stas  in  combinations 

nn  5io  kx*i.ncgnb 

ITPsNSTACfKX,  IPASS) 

LX»0 

on  520  L«t,ITP 
ITEHP«IST*C<KX,L) 

DO  510  Ksl.NSTA 
IF(I8TA(K).NE.ITEH|»)G0  to  510 
LX»LX*1 

RSTACtK*.LX,IPAS8)»K 
GO  TO  520 
51C  CONTINUE 
52C  CONTINUE 

C  REDUCE  STATIONS  TP  THOSE  IDENTIFIABLE 

N»TAC(K*.IRA39)»LX 
53C  CONTINUE 

c  identify  stations  in  tanoeh 

5*C  IF(NTNOM.LE.O)  GO  TO  400 
00  390  Ixs1»NTN0h 
DO  350  K»1,N9TA 

IFCISTAtK)  .Eil.ISTN(LXl)  GO  TO  340 
S«C  CONTINUE 
3*C  ISTNCLXInK 

NSUM(K.IPAS$)=NSMX(LX) 

ITHP«N3HX(LX) 
on  9«o  t*l , ir«P 
DO  3T0  KX9 1 »  nST  A 

KCTSTATk*)  .fo.I3TT(LX,U)  go  TO  380 
5T  C  CONTINUE 
5*C  IATfK.L.lPASslsXX 
59f.  CONTINUE 

C  IDENTIFY  PAIRS  CF  STATIONS  FOR  C0N3ISTENCT  TESTS 

40C  IF(nCSTY.LE.OJ  Go  to  43C 

■jr>  j«0  L«1.NC3JV 
DO  N 30  Kal.NsTA 

IPCISTAfKl  .EU.ISTXtLJl  GO  TO  610 
tF(ISTACK).Eo.ISTYCU)  GO  TO  620 
'.0  TO  630 
6 1 C  T9TX(L)*X 
Gn  TO  b30 
62C  I*TV(L)»K 
b3c  Continue 
6"C  continue 

65C  ITMP3NSTATNCDMB 

CD  *****  •  MAX  AND  MIN  RECORDED  VOLUMES  **•*•*••** 
C  INITIATE  8IJMS 

JM  790  KNNSTXXiITMP 
4V0DTK)»0. 

N  jaO 

DO  660  Is! ( 15 
66C  S"OfI,K)*-T 

Do  670  Tal6>  30 
bTC  SMflfI,K)«T 
TMPnO. 

TMUt*0. 

MSI 

NsO 

on  t?p  j»!,ny98 
DO  770  Isl,l2 
ISMTl 
NSN  ♦  l 

IF (X.LE.NSTA)GG  TO  TOC 
C  COMPUTE  COMBINED  710*8 

xyax-NSTA 

ItpsnST*C(KY, IPASS) 

OfM.KjsO. 

DO  690  Lsl.lTP 
ITEmP»kSTAcTk»,I.,  IPASS) 
c  CUMMTNeu  FLO*  MISSING 

IPr(!fH,IT£MP).tO.«l,,OR.Of“,K),EG.»l.)  GO  TO  680 

C  fM,K)«0(H,N)  »r.  (*',ITEmP)«CSTAC(Kx,l»IF*88) 


1213 

1219 


1217 

1218 

1219 

1220 
1221 
1222 

1224 

1225 
122* 
1227 

1229 

1230 

1231 

1232 

1233 

1234 
1233 
1236 


1238 

1239 

1240 
124] 
1242 


1244 

1245 

1246 

1247 
1246 
12*9 

1250 

1231 

1232 

1253 

1254 

1251 
1256 
1237 

1258 

1259 

1260 
1261 
1262 
1263 
126* 
1263 
1266 

1267 

1268 

1269 

1270 

1271 

1272 

1273 
127* 
1273 
1276 
127  r 


1279 

1260 

1262 

1283 


EXHIBIT  6 


GO  TO  690 

1285 

b*C 

r)(M.K)a«i, 

1266 

6“»C 

CONTINUE 

1287 

TOO 

IF(0(N,K).NlE.-l.)  GO  to  710 

12*8 

C 

staht  new  accumulations  when  flow  hissing 

1289 

N9Q 

1290 

THPaO. 

1291 

TMPAaO. 

1292 

GO  TO  770 

1293 

7 1 0 

TEMP»0(M,K) 

1299 

C 

1. MONTH  flows 

129S 

IF(SHOtl.K) .LT.TE«P)SMa(l#K)»TfHP 

1296 

IF(8H9tI*15,K) .GT.TFMPJSHfl ( J *1 5 , K ) * T£HP 

1297 

IF(5P0(n.K)  .IT.TEoPjSMCdl^iOBTEPP 

t298 

tP(SH0f?8#K)  ,r.T.TEMP)SMi;{2»,t<)aTeHP 

1299 

C 

6-honth  Flows 

1300 

TmPsTMP+TEMP 

1301 

rHp*aTHPA»TE«P 

1302 

IF(n.*)760,T30,T20 

1303 

72C 

TwPsTMP-Q (M-6,K) 

1304 

73C 

IP(TWP.LT.3M0(29,K) )8R0f29,X)aTMP 

1305 

IF(THP.GT.8Hl)(ia,K))St'  Kl4,K)aTHP 

1306 

C 

54-mONTH  FLOWS 

1307 

IF{w-9a)760.750.7«0 

1308 

7«C 

TmP43TnPA-0(m-54,K) 

1309 

75C 

rP(TPPA.LT.8r>0(30,K))SM0C3C,lOaTHPA 

1310 

IF(TMPA.GT.Jn9(15,K))SK(3(15,K)bTMPA 

1311 

C 

A  V  E  F  A  C  E  PLOW 

1312 

76C 

AV60(K)a»VCO(K)PTfMP 

1313 

N0«K*34l 

1314 

7  7  C 

CONTINUE 

1315 

7  nC 

CONTINUE 

1316 

TEHBaNO 

1317 

AVCO(K)aAVG0(W)/TE«P 

1318 

79J 

CONTINUE 

1319 

«PITEfb,300) 

1320 

HOC 

P0««aT(/J9h  MAXIMIUI  vcllpcs  pf  pecqboed  flows) 

1321 

WPITE  f6»8  10) fH0(I)»  Hi , IS) 

1322 

81  C 

rnOMAT  (5«  3TA,IZI7,33H  1-NO  6**0  54.H0  AV  HO) 

1323000 

I TNP«N3t  4  +  NCOHU 

1324 

in  A 30  KeNjTyT, ITHp 

1 32S 

IT£MpaAVr,0(i<)*.5 

1326 

JO  820  lal , 13 

1327 

02C 

I0(T)asHi)tI*K)T.3 

1328 

0  *C 

'<eTTE(6,fi«0l  I5T*(K),  CIO < I) , I»l , 15) , ITEHP 

1329 

8«C 

POPNATflx, t ".1217,218. 19, IS) 

1330 

iPITE (6/850) 

1331 

05C 

FnewAT(/iftH  minimum  vPLl'PES) 

1332 

vuITF(6,81C) CHCCI),I=1,I2) 

1333 

rn  *70  WsnsTyT, ITHP 

1334 

Do  880  1st ,  1  5 

1335 

ro(T)BSMf;(IM5,K)f.3 

1336 

07C 

■■OJTE(6,8401  T3TA  f  K  J  ,  (iQfl)  ,I»1 , 15) 

1337 

Cf  *  * 

•  •  1338 

<*ITE  (6,830) 

1339 

FOtMAT  (/21H  FseoilE’iCV  STATISTICS) 

1390 

PITE (6,890) (M0(I) , 1*1 , J2) 

134J 

«ic 

FORMAT  (/1«H  3  T  A  JTEW,I7,J1I8) 

13*2000 

c 

M[5?(NG  FLOW  PPECfOING  FIRST  NECORO  MONTH 

1343 

0 0  ooo  KBUHTyT, njta 

1344 

90C 

•J(l,H)*T 

1345 

I«CON*0 

1346 

Irfhp  a  N1TA 

1346100 

no  11*0  w*1,1TH8P 

1347000 

\r  (ITewp.CT.NJTA)  GO  TC  1180 

1347100 

K(3.LE.NST«)  GO  TO  942 

1347-2 

91C 

no  9?0  1*1, 12 

1348 

TF  WPBNl OG ( | r  W ) 

1349 

no (r,«) *oo  it, *) *,oi /tekf 

1350 

IF (00(1,8)  .LT..I)  00{I,K)a.l 

1351 

IFn-INfl»K>.LT,o.)  OC(f,w)aoO(Z,K)-OHIN(I,K) 

1352 

92C 

CONTINUE 

1353 

r>  o  o 


N»0 

93C  fin  Q«o  7al,12 
4V(I,K)»0. 

sr»(r,Kj*o. 

S*Eh(I,K)BO, 

TMPaig 

XT*CR(lJ  a(nO(I.K)+QMXN  C I  •  K  )  )/(l6.«TMP) 

9«0  CONTINUE 
942  Hal 

On  970  Jal ,NYR3 
Of)  960  lol.ts 
nbpa  i 

IF(C(M,K) .EQ.-l.l  GO  TO  930 
C  R£PU*CE  PLOW  ARRAY  WITH  LOG  ARRAY 

TFMPaAL(lG(0(H»K)+00(I,K))/2.  J023831 
tJfn.KJaTEMP 

IF(K.LE.N'9TX)  GO  TO  9*0 
c  SUM,  SQUARES,  ANO  CUBES 

AV(T,K)BAV(I,K)*TEMP 
Sn(  t,H)ssnn,K5*TErtP«TEnp 
3K6-i(I,K)*SKEW(:,K)*TE»iP*TEMP*TF.HP 
Gn  TO  960 

c  HISSING  FLOWS  ECUATFD  TO  T 

95C  ifM.KJaT 
I9C0N»1 
96C  CONTINUE 
9To  continue 

ipck.le.nstki  go  to  liao 
I  nncuO 

nn  iono  i=i. ta 
TE“P*nl0G(I»*O 
IP(TEHP.LT.3.)Gn  TO  1120 
T‘'Ps»V  (t,K) 

AV(T,K)ST*1P/TEHP 
IF(5O(I,k).LE.0.)GO  TO  980 
TMPAsSOd.O 

S0(T,R)*(St>{t,K)-4V{I,K)*TM*)/(TEMP-l») 

IF (80 ( I » X 1  .LE .0  .  )  GO  TO  9P0 

sn(r,K)»SO(I,H)»».S 

IF  («n  (I  ,K  1  ,LT.. 00(15)  GO  TO  390 

.SYFNd  pK)r(TFMF«TEMP*.«XFR  (I,K)-3.«TEMP*TnP»TMPA*2,«THP»TrtP*THP) 
1/  fTF«P*(TEi-P-l  .!«  (TEMP-2.  )*Sn(I,K)*«3) 

IF(SXF.h(I,k)  .LT.  (*.J).09.3nEw(i,K).GT..I)  1*00*1 
tF(1XEw(T,K>.GT.3.)  3XEk{T,XJ«3. 

IP  (59 EnM, A)  .LT.-1.)  SKEn(I,K)s-J. 

GO  TO  1000 
9<1C  TO(T,K)*0, 

95C  SkEW(T,K)«0. 
lone  cn'Tiw.je 

N  3  N  *  1 

IP (m.GT.I)GO  TC  1060 

<BlTP(6,lOlO)ISTA(K),(A»(T,K),I*t,ia) 

1  0  1  C  FOR-AT  (/lY,t'j,8n  Ml' AN  ,  12F8.3) 

«»TTE(6. 1020J  (SO(I,K) .  I»l#U) 

1  02C  FTN-AT  (  T  T  i  T-lS  TO  GEV,J2F,..3l 

<RITE  f6,»<*30)  (S9£p(T,a),I»I,|2J 
103C  PU-'HAT  (tOX,|MSX£N, 12*8.3) 

-=ITE(6,10«OT  (fin  ( I,  K),  l  =  i,  12) 
mne  F 'i u* h A T  (4X,6hInCm“T,F7,2,  1  1*8,2) 

■OTTf (6, 1050 T (NLOG(I,K), 1*1,12) 

10SC  POHMAT  (9K,5rtVf.ANS,l2IS) 

IO*C  IP(N.GF.IU)  1,0  TO  1180 

iP<tMnc.Le.n)  i;n  tq  lieo 

the  FCLlOwInG  BCUTINE  HILL  AOJlST  THE  INCREMENT  TC 
TRv  TO  CSTAIN  TfBC  SkEh 

Change  the  *fLLP*TN6  STAT  TO  I3*Z«1  To  ACTIVATE 

1.1X7*; 

If(TIKZ.LE.O)  go  to  1180 
ITPa-ll 

on  tuo  I»l.i2 

M  a  t  T  P  ♦  I 

On  10*0  J»l»'jTnS 


1 33* 
1355 
1336 
1 33T 
1358 
1339 

1360 

1361 

1362 

1363 

1364 
1369 

1366 

1367 

1368 

1369 
1369-2 

1370 

1371 

1372 

1373 

1374 

1375 

1376 

1377 

1378 

1379 
1379-2 

1380 
1391 

1382 

1383 

1384 
1393 

1386 

1387 
13*8 
1388-2* 
1399 

1390 

1391 

1392 

1393 


1394 

1393 

1396 

1397 

1398 

1399 

1400 
1401000 
1402 
1403000 
1404 
1405000 
1406 
1407000 
1408 
1409000 

1410 

1411 

1412 
1«13 
1414 

1413 

1416 

1417 

1418 
1*19 
1420 


EXHIBIT  6 


•6* 


Man»(A  1421 

IF  f<)  (M,K)  ,EO,T)  GO  TO  1070  1 422 

T"?aflfM,K)  1423 

'Hh.K)bIO,*«THP  *001  J,K)  1424 

GP  TO  10*0  1*23 

lore  1426 

io«o  coKTtNtie  1*27 

TFxn*SKewf I.K)  .  1*28 

IF(TE«P.flT.(-.l).*NO.TEHP.I.T.,i)  00  TO  1110  142? 

IF(Tf-Np)  1090.1110. UOO  1930 

!09c  n')(l,<)aOO<I,Kl*2.  1*31 

GO  TO  1110  1932 

HOC  O'm,K)*00{I,K)<.xtNCnm  1933 

111C  C  0*>  T  JNUE  1939 

GO  TO  930  1933 

c  *••«*«  OECFT6  3TATICNS  HIT*  (.ESS  THAN  3  YEARS  02  0ATA»  *  *  *  •  *  1436 
1120  <UJT£(6. 1130)137*00  t437 

1  1 3c  FphiaT  ( /  9  H  3T4.16.28H  DELETED,  1 NSUFF  ICIENT  0*T*)  1438000 

‘j4T*tfM)r*«!  1439 

h3T**S9ST*tI  1440 

N4T»xai'j9T**N3T*  1491 

1F(K.GT.NSTA)GC  TO  1180  1442 

C  REDUCE  S"HSC*IPT3  OF  SUBSEQUENT  37*77 ONj  1443 

on  1 1 7o  KXaK.NST*  1444 

IST*{Kx)sJSTAC-<i(*n  1443 

Hal  1446 

On  1 1  So  J’l.NYHS  144, 

On  H*0  I»l.l2  1448 

HeK*l  1449 

H'lC  Off.  KX)*3t>'»  K  X  ♦  1  J  1430 

HSr  CONTINUE  1431 

on  mo  iai,i2  1452 

7>'IH(I,KX)30MIN(I,KY*1)  1453 

M  OG ( T , 4 X ) SNLOG ( t » K X* 1 )  •  1494 

I  1  6C  OOfT, 4X)s0fl(J, KX*1)  1453 

tHC  C'STlNuE  1456 

no  rn  <no  1457 

1 1 0  C  CnTTiNDE  1438 

I T  4  *J  s  a  0  1439 

IF(TpCnH.LE.O)  GO  TO  1370  1460 

Cf.  .  *  *  .  «  AOJUflTHEST  OE  FRECUENCY  STATISTICS  TO  LONG  TERM  *  *  *  *  1461 

OP  tl?o  I ■ 1 . 1 2  1462 

OP  I  1 90  KSI.hST*  1463 

On  1 190  Lai, H STAX  1464 

NC'«4(T.4,L)=n  1463 

Si'HA(I,K,t)ao.  1466 

Si  HHCI.K.Ljao.  1467 

50A f I,K,L)30,  1468 

SO'I  f  T,X,L)oO.  (469 

XPAO  (T.K.Uao,  1470 

»T*n,K,L)S.4.  1471 

t 1 9C  PpstjnjE  (472 

DP  1220  Kcl.HSTA  1473 

xvat(»l  1474 

“=1  1473 

OP  12»0  J*1,MV6S  1476 

OP  1210  iai.12  1477 

la**|  1478 

ThmocC(M,K)  1479 

JF (TFMp.EJ.T)  GO  TO  1216  1460 

00  1200  L***,S3TAX  1481 

L»oL-nsT*  1462 

IFO-X.LT.l)  thp»;(m,l)  1483 

IF(IX.GT.O)  TMPan(H-l,(.X)  1484 

IPfTMO.EO.T)  GO  TO  1200  1483 

Nr*0(t,K,L)aMC4«(I,«,L)*i  I486 

3"*-»  (1,4,1. )’5U» MM, U«TEap  1487 

Sl,VH  ( I » 4  ,L) a;H)H4  (l, 4, 1)476?  1488 

30*  ft.4,l)aSf)*(I.4,L)*TEKPaTeHp  1469 

SOH  rr,K,L)aS04(I,4,L)*TPp*TMp  1490 

XPAH (1, 4, L) 8XBAB a, K,L )♦)£*«• Th6  1491 

IF(L.GT.nSTA)  GO  TO  1200  1442 
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*CA4{l,L.l<)3NCA8(I,Kfl) 

S'lf  *  { \  ,  L  ,  « 1 33UHH  a  , « ,  l ) 

Sll“«(I,L.K)'sattH»CX,x,L) 

3"»fI.L.K)aST»(I,K,L) 
arjn(I,L.iO=80A(J,«>l.) 
tPARd.L.KjaxPAad.X.U 
1210  CONTINUE 
J21C  Continue 
1220  Continue 
I »PC»0 

»o  t  260  K3 1  *  N8  TA 

*V3K*t 

Hi  1260  I«l«l2 

**( I.*,K)31. 

00  1250  LSK*.N9TAX 
IF(RCAlHl.K.L)  Ac. 2)  CO  TO  1230 
TF"RsKCAR(I,x.L) 

TP!jS<>-I(T,K,I.) 

T«P»e9uHA(I.K,l.> 

Tmcps  (T><p-Tf'pA**2/TE:iP)  /TEHP 

IP(T*.pp.LT,0.)  TltPPaO, 

81A(I»K»L)3TmPP«*,S 

T  -rnaSUHR  (  I ,  K  ,  L  ) 

T^^f>3(TP-Tii?s*«2/TEHpi/TEMP 

i P ( TUPP ,l  T , 0 , )  TKPP30, 

30u  ( r ,  K  , !_)  3TnPp*«  .5 

f“Pa(ff!P-TMPA*«2/TEMPl*<TP«TMPB*A2/T£hPl 
IP rTHP.Lp.a.)  GO  TO  1230 
T^«MTPAnft,K.U)-TMPA*TPpn/TEHp 
r  l  . 

!P  tTMO».i_T.0.)  TMP0o,TMP9 
TiPA3TMPA»T'1pA/r!'P 

TmPAsI .,( 1 .»TMPA)«(TEPP.1.)/(TEHP.2,) 
ip  (Tmpa.lT.O.)  T'<pasO. 

•'A(/»K|L)3THPB*THPA**,5 

(TPBl 
L  431. 

L  y  •>  L  “  N  3  T  4 

tP  (l-.LE.N9TA)  GO  Tfj  1235 
T  rp  =  t - 1 

fp  (T  rP.LT.O  I  T  PO  1  2 

L.'.-LX 

1235  IP(10(I,lO  AT.  .0001  .  OR.  SO  ( ITP ,  l  A )  ,(.  T  ,  .  0001)  GO  TO  1230 
'■i  T G  12«0 
12!r  »(TA,L)=o. 

12“:  ip  tt.  .gt.nst*)  r; n  to  i?‘io 
oOA(j,i.,x)a9  j9(I,k,L) 

L#K)sS-)A(t, y, D 
"(U.H^Ad.K.U 
IP^C  CinTJAiuE 
12NC  O'"  tTNijE 

>  12*1  *=1,nSTA 

op  I2?n  t«i,i2 

rp-'PsNLCGtl.K) 

I.  r  =  i 

•>"  1  2 T 0  L  s  1  *  n9 T  A 

Ip  (L  ,f  O.k  CI,K  ,Ll  AE.-4  .)  GO  TO  1270 

[p  (nlOGMA)  AS..M-0G{!,K)1  SO  TO  1  2  T  0 
T"P.»3NCAB(I,KiL) 

1'.P-isM.f!r.(t*L) 

T.ir.T-o»/f  1  .-(T-.P8-THPA)*K»(I,K,1.)**?/TMPB) 

IP(TP.LE.T£HP)  r.n  TO  1270 

L  y*L 

rpupaTp 

TmpPoThPA 

1 2 to  continue 

IP(I  y.UE.O)  r,0  TO  1280 

IP  (SHA  (T,K. lx)  0001. on. 80B(|, K« IX)  0001)  GO  TO  12*0 

1'iOCel 

T"P--3(3A(i,x,L*)/3nfl(i,KAX) 

T"P4S8||«An,X,L*)/THPP 


1093 

1494 

1443 

1496 

1497 
1496 
1499 
1900 

1501 

1502 

1503 

1504 

1505 

1506 

1507 
1506 

1509 

1510 

1511 

1512 

1513 

1514 

1515 

1516 

1517 

1518 

1519 

1520 

1521 

1522 

1523 
1529 

1525 

1526 

1527 

1528 


1529 

1530 

1531 

1532 

1533 

1534 

1535 

1536 

1537 

1538 

1539 

1540 
15*1 
1542 
1593 

1544 

1545 

1546 

1547 

1548 

1549 

1550 

1551 

1552 

1553 

1554 

1555 
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TMr««SuM»n,K,lX)/THI>* 

*■'  f  '  .H)eTHp*+(AV(T,l.*}-THph),»A(ItK>UX)»TMP 
3T(T  .K)»SOA(x,K,(.X)*(SO(  t»LX)»3r)«(I#K.LX))*MA(X,n,i>*j»,2liTHP 
12*0  r.^sTjuuc 

c  adjust  standard  oevtations  for  consistency 

IF  (NCSTV.LE.O)  GO  TO  1390 
C  TRAINEES  FROM  1QH 

1 2«*C  on  133o  UXsl.NCSTT 
K=ISTX(UX) 

L*ISTV(UX) 
un  1321  131*12 
TFHns(AV(I,K)-AVf I.L))/1. 

IF f AV(I,K).5T.AV(IfLJ)  GC  TO  1300 
TFPFxrEHB*SP{I,l<) 
if  (Stm, L).i.T,TEXP)  GC  TO  1310 
Tf*P=SOf t»K)#2.-TEMP 
IF  ( no  ( I  * L )  “  temp)  1320,1320,1310 
1300  TFf.paTENp*30(I,K) 

I F  t  FCi  ( 1 , 1. )  ,0 T  •  TEmP)  GC  TO  1 3 1 0 
tfkosSOT  I,K)  »2.  -  TEMP 
IF(sn(I,L).GE.TEFP)  GC  TO  1320 
131C  51(T.U*T€mP 

1 3?c  r.-sriNuE 
133C  Cl'.TIN'uE 

IFCir«v)S.GT.O)  go  TO  2820 
130c  IF (TNOC.LF.n.AKD.NCSTT.LE.O)  GO  TO  1370 
',0  T  TF  (6 , 1  350) 

13?c  FREQUENCY  STATISTICS  AFTER  ADJUSTMENTS  1 

:HITE(6,8R0) CMC(l).t=t»12) 
t'l  1361  K3l*NST» 

'<olTE(6,l0l0)ISTA(Kl,(AV(I,K),jalf12) 

xlITf  fs,112H  (son, K), 1*1, 12) 

x'0  I  TE  (6.1030)  (  SKE*  (I,  K),  1*1,12) 

•xRITEls.lOOO)  CP0(I,K),lai,i2) 

1  36C  CONTINUE 

CG  *  •  »  «  *  •  TPA  >|SF  ORN  TO  STANDARDIZED  VAR{ATCS  »***«•**«• 
13rC  DO  1«20  Kol.oSTA 

N*l 

Ofl  1410  J*1*NYRS 
10  1400  1*1  * l2 

■13N  +  1 

3R(M,K)S8LAN< 

IF  (n(N.x)  .E(3.T)GQ  TO  14C0 
I  f  ( n  n  1 1 ,  K )  ,Efj,0,)G0  TO  1390 
TfM,K)s(0()|,K)-AvtX,K))/SO(I,K) 

C  BEAPNCN  Type  III  TRANSFORM 

IF(SKEw(i,K) .EG.a.IGO  TO  l«00 
TFYP3,S«aKEM(I,K)«o(M,K)+l, 

Twp*l  , 

IFCTENP.Gt.O.JGO  TO  1380 

TP  >-Pa. tFnP 
Tm(-o-1  , 

1  3PC  Cf".’0=6.«(THP»Te>'P**n./3.)»l.)/aKEW(l,K)Y3Ken(I,K>/S, 

On  T1  1400 
13RC  3  f ,  K  )  *0  , 
me  Cn-.TINUE 
1 4 1 C  CONTINUE 
112C  CONTINUE 

C  •****•  C0"P1ITE  SUNS  OF  SGUARES  AND  CROSS  PRODUCTS  ****** 
UP  1«S0  N*l,»i3TA 
In  |<i4c  1*1,  12 
Tn  >030  l *  1 , m 9 T A x 
FAfT.X.U)  =(-«,) 

S"F*fI»F,L)SO. 

3||l*n  ( I  ,x  ,L)  *0, 

Si»  f I , x ,L)*0, 

3nr:  fI,K,L)aO, 

*®*P(I,K,L)»0. 

143c  NCA4(t,K,l.)«0 
4*(T.N,K)*|, 

1  4  4  C  CnKTINUE 
145C  CONTINUE 


1556 

1557 

1558 

1559 

1560 

1561 

1562 

1563 

1564 

1565 

1566 

1567 

1568 

1569 

1570 

1571 
1372 

1573 

1574 

1375 

1376 

1577 

1578 

1579 

1580 

1581 

1582 

1583 

1584 

1585 

1586 

1387 

1388 

1389 
1590 

*  1391 

1592 
1393 
1599 

1593 

1596 

1597 

1599 

1600 
1601 
1602 
1603 
1609 

1605 

1606 
160? 
1608 

1609 

1610 
1611 
1612 

1613 

1614 
>  1613 

1616 

1617 

>618 

1619 

1620 
1621 
1622 
1623 
1628 

1625 

1626 

1627 

1628 
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no  i9«o  ksj.nSY* 

KV3K+1 

yal 

nn  t4*o  Js 1  * nvh* 
no  1470  IB1 * i2 

TFMPa(5(N,K) 

IF(T£HP.£O.T)GC  TO  1470 
on  14*0  LB4*.NSTAX 

C  SUBSCRIPTS  EXCEED ING  OTA  BELAlE  TO  PRECEDING  PQKTH 

UX®L-M8T* 

lEfLX.LT.l)  TNP«0(M,U 
IP(t-X.GT.O)  T«P30(M-1,LXJ 
IF(THP.Ed.T)Gi)  Tn  1«6C 

c  CmNT  *N0  USE  CM.T  PECOflnEO  PAIRS 

'5uyA(i,x,l)=3'J««  ti,K#i.)*TE**p 

3H‘F\(i,K,L)3<iiiyP(IlK,l.)*TPP 

3n*  (I .K.LJaSOA  (  T  ,4  »l )  ♦TE1,P*T£MP 

3""  (1,4,1)5300  Cl,K,0*T*P»T*P 

Xf>»n(I,i<,L)»xP*cl(I,K,U*TEPP*THP 

IF (L.GT,N9TA)  GO  TO  l«6C 
*  c.Miei,u,x)5NC*»n,n,o 
3iiy*(I,L,K)s3UyM(f,i<,L) 
suy-Ml,t,K)3siiM*(i,K,L) 

90*  n,L,K)530i1  (1,4,1) 

Son  (I,l,K)=30A  ( I ,  K  ,  L ) 

X»«0(I,t,K)SxP**(I,K,l.) 

14*0  ComisuE 

IflTC  CONTINUE 
1(1*0  cn>. T T nijE 


c  *  ...»  4  COMPOTE  CUHWEL  *  T ION  COEFFICIENTS  «»*•*•*•*» 

on  i 5?o  1=1,12 

Of  1920  L=KX,NSTAX 
LX=L-N8TA 

C  ELIMINATE  PAIRS  WITH  LESS  ThAN  3  T«S  0*TA 

rE(MCAn(I,N,Ll  .LE.2)  GO  TO  1910 
TFwPSNCAH (1.4.1) 

TmP=  (  9GA  ( T,  k,  L )  *r,UN  A  ( I,  k,L)*SUMA(  I,  K,L)/ TEMP)*  (SOB  ( I,  K,L)*SUMB 
1  (I,<,L)»S0NP(I,4,L)/TE*.P) 

c  eliminate  pains  with  zero  variance  product 

If CTPP.LE.O.)  CO  to  1500 
t  “  p  n  » i , 

T''PASXPA!3  f  I  ,K,L)  -3U"A(  I  ,K,L)  «3uNB  (I,4,L)/T£NP 
c  RETAIN  ALGEBRAIC  SION 


If  (TMPA.lT.O.) TMr«a»THPP 
T«PA3THPA«THPAZT|1P 

(I  ,K,L)bTNp3*TmPA«*,9 
I  TP*  I 
L*=L 

TE(I  .IE.NRTA)  f,n  TO  HOC 
:  r  ;•  a  t  -  1 

If(tTP.LT.I)  ITPC12 

1«4C  ('<»!<)(I,K).LT.,,'P01.0S.!0(ITP,CA}.LT..0001J  H*  (1, 4.0*0, 

«o  TO  1510 
>500  RA(T, 4,0=0. 

!5!C  If  (l.  .OT.MTA)  GO  TO  U50 
-A(r,L,K)aaA(r,4,o 
1520  Cm  TIN'iE 
151c  C  r  N  T  I  Nljt 
1  5 4 C  OnKT  TNIJF 

on  m  2170 
1'jSC  ,P1TF  (6,19*0) 

1 5*0  EnRMAT(/l8H  DATA  OUT  CF  040EB) 

f,r  to  i*o 

Cw  *****  *  BE*0  COMRELATirs  COEFFICIENTS  «••«•*•«• 
INTO  DO  1*10  5B 1 »  N3T  A 

If  (K.EO.I)Gn  TO  1600 
I TPa4*l 

OP  15*0  L»1 .IIP 

c  CUWHEnT  MONTH  ccppelation 


1629 

1630 
1*31 

1632 

1633 

1634 
1633 

1636 

1637 

1638 

1639 

1640 

1641 

1642 

1643 

1644 
1643 

1646 

1647 
164B 
1649 

1690 

1691 

1652 
1633 
1694 

1653 

1656 

1657 

1658 

1659 

1660 
1661 
1662 

1663 

1664 

1665 

1666 

1667 

1668 
1660 

1670 

1671 

1672 

1673 

1674 

1675 

1676 

1677 

1678 

1679 

1 680 
16*1 

1642 

1643 
1684 

1645 
16*6 

1687 

1688 

1689 

1690 

1691 

1692 

1693 

1694 

1695 

1696 

1697 
1696 

1699 

1700 
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•*  CARD  L  **  trot 

PFAOC5.To)IT>«P,  ItEMP,  (RA(T,K.L),I«1,  12)  1702 

‘i*v(K,L)s«A(i,K,L)  1703 

IP  (TGNRL.FQ.l)  1 3  T  A  CO  *  ITfP  17C« 

IF(ITMP.NE.l!iTA(K))GO  TP  1 550  1705 

IP(mMP.M£.iSTA(L))GC  TO  1550  1706 

On  1540  1*1,12  1707 

e  RA(T,L.K)«RA(I,K,L)  *706 

C  CONTINUE  1709 

PRECEOING  MONTH  CORRELATION  1710 

C  LTSN'ST**  1711 

Ip  f IGNRl «£0 . t )  lX»N3Ta»k  1712 

LAsnsT  At  1713 

IP  flGNPL.eO.l)  LA*L*  1714 

M  1610  L8L*,LA  1713 

I T  p  *1  «HST*  '716 

•*  CABO  K  OR  H  ••  1717 

HE *0(3,70) IThP, l TEMP, (RA(I,k,L) .Is  1.12)  1718 

IF  (K.F-O.lUSTAtKJsiTiiP  1713 

IP{R.E0.DI9tA(ITP)sITEPP  1720 

IP  (IGNBL.EO.t)  RAV(K,K)cBA(1,k,L>  1721 

IPtITMP.NE.I",TA(K))r,o  TC  1530  1722 

l P  C I TEmP •  N£ , t  S  T  4 ( I  TP) ) GO  TP  1550  1723 

3  CONTINUE  1724 

on  1620  I»l»l2  1725 

3  'H(T,K,K)al,  1726 

3  continue  1727 

*  .  *  *  •  READ  PHEOuENpY  STATISTICS  *«***•*•*•»«*•**  1726 

On  lbflO  K=1#nST*  1723 

,  **  CABO  N  [P  0  »•  1730 

'<  PA  0(5, 80)  ITP.(AV(I,K),I*1,12J  1731 

rp(tTP.NE.I.4TA(K))G0  70  1550  1732 

GENpRALl7tO  STATISTICS  On  ONE  C*RO  PER  STATION  ITJ3 

Avxy  (K)s*V( 1 ,K)  .  1734 

a  v  n  n (k)bav(2,K)  1735 

snAv  (AJSAV  (3,K)  1736 

I  rPPnAVf'WK)*,!  |  737  * 

I"Y  f5)SITMP-f»0(12)  1736  • 

[TuPaAv(5>X)6*l  1738-2* 

I"N  (;<)eTTMP-MI(12)  1739  * 

ip  (IHX(K)  ,LT.l)  I»«*(K)»IP»(K)*12  1740 

IP  ( I  NO  t K ) ,UT • 1 )  lMN(K)SIPM(K)Tia  17«1 

IP(HJN«L.E0.l3Sn  TO  1 6 *10  1742 

•  •  CABO  P  *•  |7«3 

BFAP(5,S0)ITM,(.9O(I,K),I*l,12)  17*4 

IP(TTP.NE.ISTA(X))00  TO  1550  17*3 

•*  C»SO  0  ••  1746 

REAn(j,#n)jTp, (SKEwa »*),!«! ,12)  17*7 

IP (TTP,NE, JST A (K) ) 00  70  1350  IT*8 

•*  CANO  B  •*  |7*9 

!iPAO(3,90)lTP,(OOU»A)#t»l*12)  1750 

IP(ITP.NE.ISTA(:<))G0  to  1550  1731 

3  CONTINUE  1732 

*  .  «  *  *  ESTIMATE  MISSING  CORRELATION  COEFFICIENTS  »•***•••  1753 

)  IP ( IGNBL .EC. 1 1  GO  Tn  3020  1739 

IF  ( n jt A ,i.K ,  l )  GO  TO  231'i  1755 

00  1720  I3 1 r l 2  1756 

IPSI-1  1757 

IPITf.tT. 1)10=12  1756 

!)0  17)0  K  =  1.nSTA  1759 

I T  P  s*K  ♦  1  1760 

On  1700  LaiTP,NSTAx  1761 

L  AnO  k  CORRELATION  P0SSJ8LT  MJ33IN5  1762 

IFtRAfX.K»l>.6E.C-l.)>  GO  TO  1700  1763 

RmatbI,  176* 

B  m 1  wo* 1 ,  1765 

ux  SEARCHES  all  fltlATEO  CORRELATIONS  EXCEPT  FOLLC*INS  ATM  1766 
On  16*0  L**1.NSTAT  1767 

IF  (1.  v.EO,R)Gn  TO  1690  1768 

IFCL.EO.LX)®'1  Tn  ,6,°  1763 

TP»Dr(JA(I,K,L<)  1770 

1F(L,LE.N1T*)GC  TP  1660  1771 
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IF  (LX.LE.NSTaK.CI  TO  1670 

c  8otm  l  and  lx  f/efirfsesr  p«eccoi*8  hont*t 

ITnpaL-NaTA 
tTEMP«l,X«NSTA 
T»‘P»RA  {II**  ITmP*  I  TEMP) 
on  TO  1680 

C  L  0EPFE88HT3  CURRENT  month 

|*A0  T"PaR*(I,L,Lx) 

Gn  TO  i860 

c  LX  ANO  NOT  U  REPBMJNT8  CURRENT  MONTH 

ifcTO  T*'PsBA  ( I #LX i L ) 

1 6*C  IF  f TPP*tEHP.LT,«2.0)  GC  Tr  1600 

T-PAaffl  ,-tEmF*TEnP)  *(1  ,*TMP«TpP)) 

IFfTMPA.LT.O.)TNpAsO. 

TmPAsTmPa**.5 

T>'?RsTHP*T£Mp*TMPA 
IF(THP«.lT,Bm*  V)P><AX»ThP6 
rppMaTMPS-TMBA-TMPA 
IF  rTMPfl.GT.PftTOKMlNPTNPP 
16*C  CONTINUE 

c  AVERAGE  SMALLEST  pax  and  LAPSES  I  H IN  CONSISTENT  VALUE 

BA  CT,K,LJ  =  (PHAX»N«*!N)*,S 
IF (L.LE.NSTA)PA (I,L,X)»flAfI,K,t) 

170C  cn» TINUE 
1 7  l C  CnsTINIlE 
1 72C  C^rTINue 

Gt  TO  2313 

cj  *****  *  test  for  triad  consistenct  ************* 

I  7  1C  NCAo<j 
17«C  FAC*1, 

NO  asnC  A*1 

IF (NCA.LT.NSTAF12)  go  TC  1750 
.-.t?  I  TE  Cfj  *  1 8403 
On  TO  150 
1 75C  MV<so 
NCs.1 

1 760  { >  TCsQ 

On  l»TO  I»l « 12 
tP"T*l 

TF(IP.LT.niPal2 

c  K,  L ■  a kO  LX  SEARCH  ALL  BELATED  THlUS  OF  COMPEL  CCEFS 

on  <*20  K*1 1 N8T A 
tT'P=K*l 

OH  ISJO  L=ITmP, NSTAX 
(l.Etl.NSTAx)GO  TO  1*10 
L «  fI.-NSTA 
N  t  a  i*  A  f  1 1  K  a  L  ) 

I T  USL ♦ l 

On  HO0  L*  =  tT°.  "STAX 

ITF»P3LX«N3fA 

P,’SI7A(I,K,LX) 

IF  (L.LE.N5T*)RJ!5“4(I»l»tXl 

c  0 n T ■<  L  AND  L»  REPRESENT  PRECEDING  MONTH 

IF  CL.GT.NSTA)R3sB*(If*«LA,  I  TEMP) 

C  RAISE  LOWEST  COEFFICIENT  IF  INCONSISTENT 

ACts (| ,-Rl*Bj ) *o,3 
A  C  2  s  ( 1  ,-p2»P2)»*,5 
AC «s(l ,.SJ*R1)**,3 
IP  C -» 1  .ST.S2)  GC  TQ  1770 
IFtPl.GT.B3)  O.G  rn  17*0 
7>'  i  *  sP£*P 3* *C?  « *C 3 *FaC 
IFfp-TN.LT.-l.)  Rmxns.I. 

IF{M1,GE,BmI*j)  Gn  TO  1*00 
I»  ncet 

RA ;I.P,L)*»HIN 

IF  fL.LE.N3TA)  PA(I#L.K)B»HIN 
GO  TO  1«00 

17TC  IF  C’JP.GT.BJ)  GO  TO  1710 
t<Mj^am*B3-ACl«AC3«FAC 
JF (J"IN,lT.*1 .)  PMINP.l. 

IFfea.GE.BPIN)  on  TO  1*00 
INOCBl 


)772 

1773 

ITT* 

177J 

1776 

1777 

1778 
1773 
1780 
170] 

1782 

1783 

1784 
17*3 
1706 

1787 

1788 
1783 

1730 

1731 

1732 

1733 

1734 
1733 

1736 

1737 

1738 
1733 
1600 
1601 
1802 

1803 

1804 
1803 
1806 

1807 

1808 
1803 
1810 
1811 
1812 

1813 

1814 
IBIS 
1816 

1817 

1818 
1*13 
1820 
1821 
1822 
1823 
1*24 
1823 
1826 

1827 

1828 
1*29 
1*30 
1*31 

1832 

1833 

1834 
1833 

1836 

1837 
1*38 
18)9 
18*0 
18*1 
18*2 
1843 
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ha(t,k,lx)»Rmin 

1844 

IP  (LX.LE.NST*)  RA(1,LX>K)»RMI* 

1845 

t;n  to  i*oo 

1*46 

1780 

i)MINeRl«R2-ACl*AC2«FAC 

1847 

IP(9HIN.LT.-1 ,)  HMINB.l, 

1848 

IFfB3.GE.RMlN)  GO  TO  1*00 

18*9 

INOCa! 

1850 

IF  fL.GT.NSTA)  GO  TO  1790 

test 

RA(t,L,LX)eRMtN 

1852 

IF  fLX.LE.N3TA)  HA(I,LX.L)«RnXN 

1853 

on  TO  1*00 

1854 

179c 

BA  (TP.LA, ITEHP)aR8XN 

1*55 

PA(IP,IT£HP,L*)aRMIN 

1656 

tunc 

Continue 

1857 

i  a  i  c 

C  n  8 tinuE 

1858 

t  82c 

Continue 

1*59 

133C 

continue 

1860 

NCaNCM 

1861 

IF(NC.LE.N3Ta*1?)  GO  TO  1850 

1662 

'-51X6(6,18401 

1663 

18«0 

fokm«t(32h  Correlation  matrix  inconsistent) 

1864 

SO  to  150 

1*65 

10*0 

IF(INDC.EO.l)  GO  TO  1760 

1866 

r,<  •  * 

*  .  *  •  TEST  FON  OVER-ALL  CONSISTENCY  *****•«• 

ir£np»o 

1868 

f,o  TO  1«70 

1869 

106C 

ItEmpbi 

1870 

when  i temp°i ,  current  month  used  r as  all  independent  stas  i at i 

OTHERWISE,  PRCC  "TH  USED  FOR  CURRENT  AND  SUBSEQUENT  STAS  1872 
‘4  T  KOP3NST  A  1673 

:,V*  NsNINOP+t  1874 

ini  2160  1®1,12  1875 

tpst-l  1676 

leUP.LT.nin3lZ  „  •  »8?7 

CO^S  TSiJC T  COMPLETE  CORREL  MATRIX  FOR  EACH  MONTH  AND  8 TA  1878 


1  A^C 
ibrc 


I90C 
I  9lC 


t  9AC 
!9lC 


19  4.; 


t  ISC 
19*C 
I  970 
1  RAC 


199C 


2ooe 


0"  2150  K  =  1 » *!ST  A 

L  IS  sow  NUMBER,  J  IS  COLUMN  NUMBER 
.in  3020  Lsl»NSTA 
Lr*L*N3T* 

On  1930  J»t,NSTA 
J »s JaNJTA 

Tr  (i  -K)  1080, 1920,1960 
(r<!.K)  1890.1910. 1900 
9fL. J)  =  njlE(AA(I,L,JJ) 

I.TMB(L)»l 
JT'P(J)«J 
r,o  rn  1970 

IF(ITEmP)  1910,1910,1890 

-a.j)  »  ohle(na(i»l#jx)) 

LTM«(L)«L 

JTmp(J)*JX 

■ir  rn  i9T0 

IFf.r.K)  1930,  1990,1950 
S(l.J)  ■  ORLE(RA(I,J, un 
LThp(d«j 
JTPP(J)*L* 

on  rn  1970 

«  DOL£(RA(TP,L»J)) 

LTPPa)«LX 
JT“P(J)«JX 
GO  TP  1970 

U(fTffP)  1990,1990,1930 
IP  f  TTfMP)  l9i'0,  1920,1880 
of j,l)9R(L,J) 

CO' TjNyC 

iv»ai« 

SPECIAL  SUBSCRIPT  FOR  DEPENDENT  VARIABLE 
IP  (L-A)  199n, 2010,2000 
sa,*49TAA)  ■  OPLE  f  HA  (I,R  ,L)) 

J  f^u(NSTAA)»(. 

50  TO  2020 

IE  (ITEMP.CT.O)  f.O  TO  1990 


1879 

1880 
1881 
1882 

1883 

1884 
1685 
1886 
1887000 
1888 

1889 

1890 

1891 
1892000 
1893 
1899 

1895 

1896 
1697000 

1898 

1899 

1900 
1901000 

1902 

1903 

1904 

1905 

1906 

1907 

190e 

1909 

1910 

1911 
1912000 

1913 

1914 

1915 
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EXHIBIT  6 


ao 1 C  Rfl.,H8T*A)  °  oble (R*(I,K,LX)) 

JTH»(NST**)SL* 

202C  CPNTtNuE 

C  MATRIX  CONSISTENT  IF  CORREL  ODES  NOT  EXCEEO  t.O 

NaO 
NCaP 

C  KSZ3SS1B3S 

ao30  CAu.  CROUT(R) 

C  B=3E33aasS 

1F(DT»MC.LE.1.)  GQ  TO  2130 
RRITE(6,2O40)  N,I,k,OT»nc 

?o«c  format  r/36H  inconsistent  correl  matrix  adjusted, 3I4,fi2. 3) 

c  ,  kithcrah 

fac*fac-.2 

IF  (FAC  .GT«**»  1 )  GO  TO  1750 

NCOal 

nbn*| 

IF(N.DT.lO)  GO  To  150 

StlMsO , 

on  ?0A0  l ® 1 , N I NOP 
no  ?o7o  u***«nvar 
IFCL.E'J.LX)  GO  TO  2070 
TNPOsR (L.LX) 

S'IM3Si)M*THPF 
207c  CONTINUE 
POFC  CONTINUE 

TFmp  =  nINO**ANINDP 
S'inssum/TEmP 
T  E  ^OaOTWFC* l . 

IFCTFHP,r,T,.D  TEMpa.l 
TMPal ..TEMP 
00  P120  L  *  1 i NlNOP 
ITPbi.  h 

no  7U0  i_XalTP,NvAR 

RfL.LXJ  =  OBlE (ToPP«TNP  ♦  SUpotENP) 

IFax.UE.NINoF)  P(l.X,l)sR  (L.LX) 

LTVauTrlPfL) 

JtrsjthP  (LX) 

IFtt.TP.LE.NSTA)  go  to  2100 

IFdTP.Lt.NST*)  GO  TO  2090 

LTFaLTp-NSTA 

JTU3JTP-N9TA 

RA  f  TP.l.TP,  JTp)aR(L,LX) 

RMIP.JTP,LTP)=R(L,LX) 

Go  to  2 1 1 0 
?OlC  ITMPCI.TP 
LTPa.ITP 
jTPatTMP 

2 1 OC  TA(T,lTP.JTP)3A(L»LF) 

IF  (fTP.LF.N3  TAJ  NA  (I,.ITP,UTP)=R  IL.L*) 

a i » r  continue 
2  1  2C  CONTINUE 

in  to  2030 

a 1 3C  IF  (OTU  IC.GE.O.)  GO  to  2140 
-9ITE(i,,7OJI,K,0T»(!C 
OTP ”C*0, 

a t  - 3C  IF (NCfl.GT.O)  GO  TO  1740 
2 1 5C  continue 
?  1 6  C  Continue 

IF  f ITEMP.EO.O)  go  TO  I860 
if  (ITUNS.EQ.2l  Gtl  to  3100 
?  1  7C  •<«ITffs,l01 

C  •  PRINT  CORRELATION 

on  2260  I»l,l? 

IFf  TTRNS.LE.ft)  WG  ItE(6,21*0)m0(I) 

21^0  Format  I//39H  R A h  CORRELATION  COFFFICIEnIS  FOR  MONTH, 13) 

IF  (  TTRN8.GT  .0)  Hn’ITE(8,?190J  hoTH 
?19C  Fofhat  (//90h  CONSISTENT  CORRELATION  NAT6JX  FOR  MONTH, J3) 

»OITt (6.2200) (1ST  A (K) ,K«l,N8TA) 

2230  Format  ( /3X , 3HST A , 1 81 7 ) 
mpITE (6,2210) 

2210  F  Ohm  at { 20X , 1 9H  aITH  CUhpENT  MONTH) 


ITR  l ) 


1916000 

1917 

1916 

1919 

1920 

1921 

1922 

1923 

1924 
1923 
1926 
1927000 

1926*1931 

1932 

1933 

1934 
1933 

1936 

1937 
1930 

1939 

1940 

1941 

1942 

1943 

1944 

1945 

1946 

1947 
1946 

1949 

1950 

1951 

1952 
19330O0 
1954 

1953 
1936 
1957 
1950 
1939 

1960 

1961 

1962 

1963 

1964 

1965 

1966 

1967 
I960 

1969 

1970 

1971 

1972 

1973 

1974 

1975 

1976 
t  977 
1970 
1979 
I960 

•  •  1961 
1902 

1963 

I  984000 

1963 

1986000 

1967 

1980000 

1949 

1990 


•  1 4« 


EXHIBIT  6 


22?e 

22  50 
22#C 


2250 

2a*C 


C 


22TC 

?2*C 

22«C 

2300 


HO  2220  K«i,i,STA  .  ...... 

'"me <6, 2230 >  1 

enSM4T  (t*,t5,l«rT,S) 

format  (iox,36H  hitm  PhfCEOI^e  month  at  above  station) 

ITPoVST**! 

on  2250  KBJ.NSTA 

-:0  TTf.(<,,2230)lSTA(K  ),(»»{ I  ,K  ,L  )  »L«  I  TP»  KST  AX  J 
Continue 

IF(IANAL.te.O)  GO  TO  HOC 
IFfir»NS.LE.O>  or  TO  1630 
miPCHS.u.o)  nr  TO  2670  „ 

2U*'CH  ESSENTIAL  elements  op  matrix 

Cn  2300  K»l,"STA 
IP  fK.EQ.t)  GO  TO  228C 
T  tps<- t 

"o  2270  L»l.lTP  ... 

mTF(7#7<')nT*(K),  ISTA(t)»  f9A  ( I#  Ki  L)  »  I*  *  » 


no  2260  L«n9TA4,|.0TAX 
TTF1P3|_*NSTA 

*«ITKf7,70)  TST4(X), J8T*{ITEmP), 

Continue 


£RA{I<A,'-),I»J  »)8) 


GO  TO  2850 

Cl.  *****  *  RECnN.>iTtTuTE  MISSING  0*T 
2310  IFfTANAL.LE.o)  GO  TO  3 1  <*C 
IF  ftGcrN.LE.O)  r,n  TO  2*10 
':V*PSN3TA*1 


!<3l 

On  ?600  JB1.MYH3 
On  2560  !■! t 12 
lYST-1 

IFdX.LT.  1)17*12 

M  =  M+t 

DO  2500  K a 1 1 N S T A 
IF{GfM,K).NF.T)  GO  TO  21*0 

C  F0«*i  C.OHRELAT  If  N  MATRIX  FOR  EACH  MI88INS  FLO* 


K  f NOOaO 

on  2J63  LPlflStA 

L  v -I  *N3TA 

tf{o(m,l)  ,nf.t)  r,n  to  2120 

IF(0(M-l,t).r>,T)  GO  TO  2360 
\T».np«fjiNnp*i 
I Tf MPBNTNOP 
X  (► IN0P}*O(M«l ,L) 

#(MnOP,NVAH)  »  O'iuECRAdfX.l*)) 
Gn  TO  2330 
2 J2C  ntnOPbnIN 02*1 
ItEmpsmTNOP 
XI"  tnOpJ  *0  (M ,  L  ) 

*.  ?‘-TM5P.NVAR)  S  tJ?LE(»A(I»F,L)) 
233C  5  t  '■  !►  Op.t.  JM1P)  *  l.OOC 

rc  n.  .eo.nsta)  cn  to  216C 

ITF3|.*t 

pr  23*0  L*«ITP.NSTA 
Jyai  t ♦  n 6 T A 

IF (nfM.L) .EO.T)  GO  TO  2350 
1F(o(m,i. a).EO.T)  GO  TC  23*0 

l TF MPSITEMP*! 

A  ( *  t  MOP ,  I  T£MP  5  3  OBLE(RAd.L.LA)) 


GO  TO  2370 

2  3  JC  IF  (.»fM-l,LA).C5.T)  GO  TC  23«0 
ITF-p»ITf MP*l 

Pf‘ Tf  OO,  fTE"Pl  »  OHLEfRAd.l.JX)) 

GO  TO  2370 

231C  IF  (1{m,LA)  .EiJ.T)  GO  TC  2360 

WfNIMnp.iTEMP)  ■  OFLE (R* {I«L*  »LX) 1 
GO  TO  23TO 

2J6C  IFfO(M-),LA).EG.T)  60  TC  *380 
I TF  hm* IT£MP* 1 

RdTNOP.lTfMp)  a  r,0LE(RA(l7»L.lA)) 
AOO  SYMMETRICAL  ELEMENTS 


1661 

1662 

1663000 

166# 

1665000 

1666 

1967 

1668 

1966 

2000 

8001 

2002 

2003 

200# 

2009 

2006 

2007 

2008 

2009 

2010 
2011 
2012 
2013 
201# 

2015 

2016 

2017 

2018 

2019 

2020 
2021 
2022 
2023 
202# 
2023 
2026 
2027 
2020 

2029 

2030 

2031 

2032 
20J3 
203# 
2033000 
203b 

2037 

2038 

2039 
20«0000 
20#1000 
20#2 
20#3 
20## 
20#3 
20#6 
2007 
2O«0 
20«9000 
2050 
2031 
2052 
20530C0 
205# 

2055 

2056 
2057000 
2036 
2039 
2060 
206 1000 
2062 
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EXHIBIT  6 


33TC 
33*C 
3  39C 


340C 

3U  JC 
C 

342C 

C 


?l!C 

3««n 


c 

c 

3«*i0 


?u6C 
?  1  7  C 

?'!lc 


C 

35K 


Ji?0 

c 


353C 


C 

C 


C 


RfITEMP«Nl(JDP)aR(NINOP,IIEHP) 

2063 

CONTINUE 

206® 

CONTINUE 

2063 

IF(NIM)P.GY.O)  GO  TO  2400 

2066 

MNOP«t 

2067 

xin*o. 

2066 

pfi.n  ■  i.ooo 

2069000 

LXaK*NSTA 

2070 

Pfl.NVAR)  3  D8lE(RACI,K,LX)) 

2071000 

I  r E hoon I  NOP ♦ 1 

2072 

(in  ?U10  L*t » NtNDP 

2073 

UfL. rTEMPl=9(L»NvAH) 

207* 

2073 

call  crout  (pi 

2076 

:::::i:333 

2077 

ITEupaNlNOPn 

2078 

TFPP«1  . 

3079 

InDCsO 

2080 

On  34*0  l  3 1  .  n  T  n  0  p 

2081 

ThPBO  AR3  (R  (L  i  HEMP)  ) 

2082 

IF(Twp.GT.TEnP)  GO  TO  2*30 

2083 

rpM PsTmP 

3084 

T  TP»L 

2083 

IF (P (L » ITEhP) ,LE.0..AN0.atL).«T.C-t.5).AN0.B(U.LT. 

0.5)  GC 

to 

3*40  2066 

IP (9(L, iTEMPl .GE.0..AN0.8 U) .Gl . (-0 .5) , ANO .6 (L) ,LT, 

1.5)  GC 

TO 

2**0  2087 

iNOCal 

20*8 

CONTINUE 

3089 

IF CTnOC.GT ,0)  GO  TO  3«50 

3090 

IP  CDTRMC.LE.  I  .  .ANP.OTP'iC.GF.O.)  GP  TO  2310 

2091 

IF  mA  rpj*  INCONSISTENT.  OMIT  VARIABLE  *ITH 

LEAST 

2093 

cupnelation 

3093 

JT"P*NINPP.» 

2094 

IF(ITP.GT.JTmP)  ca  TO  2480 

2053 

on  ?'iTo  l  =  ITi*.  itnp 

3096 

no  ?*6o  L*si .  itemp 

3097 

PfL.LA)»»(L+l*LA) 

3098 

Vfl  laX(LTl) 

3099 

i'p  3500  L a  1 .  I T N P 

3100 

Up  3«P0  L*=lTP,NINnP 

3101 

a.LAl=HlL.L*Tl) 

3102 

Continue 

3103 

*j  T  NOpa  I TMP 

3104 

G"  TO  3930 

3105 

AOD  PANOOM  CPNPONFNT  TO  PRESERVE  VARIANCE 

2106 

I p  NP30 . 

3107 

0.1  3520  l*1.N 

2108 

TP“P  =  TEMp*hNi,En  ( I X X ) 

3109 

TPvOaTENP*PNGEN(I*X) 

3110 

compute  flow 

2111 

*1  3  f  1  ,-OTP‘iC)  *«,5 

2113 

TPPP3TEMP*AL 

3113 

on  ppjo  La 1 »  NlNOP 

311* 

TrK9aTFMp*n(L)«*(LJ 

3113 

"'“.KlaTEMP 

2116 

(.:•  tM.Kjaf 

3117 

ro=OfM,K) 

3118 

ADD  NEN  VALUE  TC  SUMS  OF  SQUARES  ANQ  CROSS 

products 

2119 

Or  3560  1»J.*J3TAX 

3120 

IF(L.EO.K)  GO  TO  3560 

2121 

subscripts  e»ceeoing  nsta  relate  to  preceoing  numh 

3122 

Lvai  -NSTA 

2123 

TFtLX.LT.lJ  TMPaoIM.Ll 

2124 

ifux.gt.0)  TMP=r.lH-i,ix) 

2123 

IF  tTNP.EO.T)  GP  TO  3560 

2126 

COUNT  AhO  USE  CM*  RECORDED  PAIRS 

2127 

•<P  *  0  c  1 ,  K  ,  L  )  =nC  A8  (  T  ,  *  .  U  ♦  1 

2126 

MlMAfI,K,L)aSUNAtI»K.lJ*TP 

2129 

5  IN'«(I,K,U«SUNB(T,K,L)*TNP 

2130 

ST  A  (I,K,l)asOA  (T.R.L)v  T  R»TP 

2131 

So»  (I.N.LjasOO  CI.R#U*7NF»THP 

2132 

»PA0(I,*,L)«<FA8CI,K,L)«TF»TMP 

2133 

if  Cl. St. NSTA)  GO  TO  2540 

213* 

•  It* 


EXHIBIT  * 


‘jcJMUji.KiBMCAin.K.u 

Sl'N*tt.L.tO«3i)>-R{I#«.L> 

?ltK«fI,|.,K)03tJI'A(t,K,L) 

JOA  (t.L.K)asnM  (t,K,L) 

SOP  (I.L.KJCSQA  <I,K,L> 

XPAB(I,L,K)axPA8<I,K,l) 

c  recompute  correlation  coefficients  to  Include  nea  oat* 

?5«0  IR(nCA8{I,k»L) .U.2)  GO  TO  2560 
TEhp»nc*B(I.K#D 

TMPr(80A(I#K,U»SUMA(I,K,L)<‘SliHA<l|K»L)/TEMP)*<a8B(I.K,L)'*3UH8 

l(I,K,L)*3UHBa,lt,L)/TEMP) 

c  eliminate  pains  «ith  zero  variance  product 

IF(THP,L£.0.)  GO  TO  2560 

THPilsl  . 

TMPAaxpAB(I,K»L)-3U'lA(I,K#U«aUNBtI*N*l.)^TEMP 

C  RETAIN  algebraic  sign 

IP  f  TMPA.LT. 0.)  THpfts-TNPS 

T*<PAsTHPA*TMPA/TpP 

RAC  7 ,K,L)  =  TNPn*TMPA»*.5 

ITPOI 

L*aL 

IPCL.LE.NflTA)  GO  TD  2550 
I  TP  3 1 • t 

IPCTTP.LT.l)  I T  Pfi 1 2 
L  A  *1  V 

2550  IPfO(I,K).LT..OO0l.0».S0nTP<LA).LT..OOOl)  RA(I,K,l)aO. 
IPCL.5T.N9T*)  GO  TU  2560 
C» CT,L,K)3RA(I,K,L) 

25*0  CONTINUE 

ITiP3NvRS«12a1 
IF(m.GE.ITmP)  GO  TO  2530 
TP‘ p=C (M.K) 

OP  ?«7o  L*l *  U9TA 
T  “PB')  (f*+  I  f  L  ) 

IFCTttP.EO.T)  GC  TO  2570 
LV3**n$TA 
I  TPs 7«l 

lF(TTP,f)T.j2)  I TP3 1 
NT  IP ( JTP,L,LX)SNCA«{ ITP.UfLX)*! 
!5i'”AfI'fP.L,Lx)=9uJACITP,L,LV)tT,»P 
GUM A(ITP.L.LV)s9u”B(lTP,L»l X)*TP 
SOA  (ITP.L.LXjsSGA  CITP,L.L*)*TMP*TPP 
S19  (ITP.L.LX)SSU9  (ITP.Lil *)*TR*TP 
7P1«  CITP.L,Lx)=XP»dCITP,l,LX)+TP«THP 
IPfRCA!»(lTP,L.L*).LE.2)  GO  TO  2570 
TPMPSMCA6 CITp,L.L«) 

Tops  f oqA  flTP.L, Lx) -SUMAC  I  TP, l ,LX) *SUM * ( I  TP , L , LX ) /TEMP)* 

1  (9i;>C  ITP,L.LX)-5tiM'»(ITP,L»C  X )  tjUMP  ( I  TP,  L»  EX) /TE*P) 

Ip CTVP.LP.O.)  GO  TO  2970 
TMPA=l  . 

TMP43XP»aCITp,L  ,LX)-3LM» (ITP,L,LX)»5UNB(ITPrL»L*)/TEMP 

IT  (T'tPA.LT.O.)  TMP9a«TMPB 

TmpAsTmOAATOPA/TmP 

t:‘crrp,L.Lx)=TMPB*T*<PA*».3 

1PC80(I,K).LT..OOCI.D«.SO(ITP,L).LT..OOOI)  RA C I  TP , L > LX ) »0 , 

’570  CP‘ TfNue 
?5*r  r.nsfTNue 
2'5qC  C  'NT  I  nijE 
?<>OC  C  pn  T  I  nuE 

26 1 C  IP  CIANAL.LE.O)  GO  TO  3100 

CM  «  »  *  •  «  •  CqnvEsT  3T»nCAS0  revIATES  to  FlQ«3  ********* 
IP  C'lOAgA. LE  ,  l )  Gu  TO  2430 
ITMPpe>iroS»l?»l 
On  2620  I Tt*P»  1 , ) 00 
IPCIOTAP.EO.mOTAP)  GO  TC  2630 
HP*n  f  TOTAP) 

262C  L0T*B3L0TAP*1 
263C  <PITE(6,!0) 

mb  lit  (6,26*0) 

?6«0  PPnMAT(33M  RECCROEO  ANO  RECONSTITUTED  FlQn») 

IF ( ‘ PAG5.GT  .  J )  "R I TE (6,2650  !P*SS 
265C  FPRMATC5H  PASS, 13) 


2135 

2136 

2137 
2136 
2139 
2l«0 
21«1 
2162 

2143 

2144 

2145 

2146 

2147 

2143 

2149 

2150 

2151 

2152 

2153 

2154 

2155 

2156 

2157 
2156 

2159 

2160 
2161 
2162 

2163 

2164 

2165 

2166 
2167 
2166 

2169 

2170 
21*1 

2172 

2173 

2174 

2175 

2176 

2177 
2176 
2179 
2160 
2161 
2162 
2163 
2184 

2165 
2186 
2167 

2166 

2189 

2190 

2191 

2192 

2193 

2144 

2195 

2196 

2197 

2198 

2199 

2200 
2201 
2202 
2203 
2206 

2205 

2206 
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Exhibit  % 


ANYBSSMVSS 
pn  2*10  KS1 . 'iSTA 

ip (k.st.kstx)  KRirr(6,ai!.(,c)  (Modi.iaid2) 

266  C  FnOMAT  C/llH  ST4  T£  A  » ,  1 2  ]  * ,  t>X  #  36  TOT  At.? 

Hat 

'in  2760  J*t#'JVBS 

ITPaO 

on  2720  1*1(12 

MS6*1 

TF«Pa<J(Nf  KJ 
TmpsSKEHM.X} 

IFfTMP.NE.O.)  TENPa((TMF«(TFMP.TNP/6,)/6iM.)**3  -l.)*2./TPP 

iFcopfM^Kt.'JF.n  go  r r  ?e«o 

IF  (tEmp.GT.2.  .Anc.SDCI  ,K)  ,r,T.  ,3)  T£MP»2.* (TEMP-2,) *.3/SQ( I (K) 

IF (TKP.LT.-.onot  .09.TFP.GT..0001)  TMP»(-2.)/THP 
IF<akEw(I,k11  ehTP.abSO.JfcttO 
2&7C  IF  (  rf-P.GT.T  iP)  TEMPaTFP 
'in  TO  2690 

26*0  IF (TEKP.lT.TmP)  TEhPbTOP 
2690  T',paTEMP«SOfT.K)»*vCX»K) 

'!  (P  t  K 1  ■  1 0 .*»  THP-UO  ( I  »K  ) 

iF(a(M,K).ur.n..ASo.oFiK(i,x).Ge.o.)  aCM»Kj»o. 

.«(I)n#(H,K) 

ITHPawSIlHlK,  I  °  A  S  3  ) 

IFf  ITMP.LE.O)  GO  TO  2?10 
TF“PaO. 

in  2700  Lat.IT“P 
I.VST4T(K,L(IP*S9) 

27  00  TF^P«T£MP*0(m.IX) 

IF  ( ')  ( M  t  K  )  .GT.TE-Pl  go  TC  2710 
r.*  I  t)*ADJ 

TF(00(m,k).NE.E)  GO  TC  2710 
tnaAOjl 

:;(“(K1*te**p 

27  1  C  I'KT)«g(N,K)  +  .5 
2720  ITF=ITP*IO(n 

;yc«|YpA»j 

■F (K.LE.NOTX)  Gn  TO  2760 
IF  (TPCHO.LF.Olun  TO  27<|0 

7TF I7#27 JO)  I3T»(K) , ITS, (toll), la|,121 
2710  FP»HAT(2j«,iaib) 

p 7 a C  ;iTE(b,2750)rsfA{>O,IrP(  tIQ(!)(QPU)d«l(l2)  dTP 
3710  Fr":nAT(lX(I«d6d8(M»ltCI7»M)dtO) 

2760  CPatinuE 

IF (UPASS.LE.  1)  GO  TO  2765 

>ntTEfIOT»P)IU*tK)((C(i',*).Pai(lTMPPl 

-,n  iapswotapm 

3765  U-  I  f  9C0N.LS.0l  GO  TO  281C 

CK  *  •  *  *  •REcnnai'TC  MEAN  A(,r  3T  A  NO A  90  DEVIATION  •  •••***•»* 
>n  2770  1*1.12 
AV(T,K)»0. 

GIF <(I,K)SO. 

27 7 C  jO  (  I  ,K  )  »0  . 

'cl 

in  3790  ,1  a  1  .  N T F  3 
in  37*0  J ■  1  >  1  2 
n  =  r-  +  I 

TFFPeAL0P(0fM.><)*00fI,X))».P3«29fl3 
aw(T(K]bav(  I,«<)*TEhP 
i/Ht  /(I.<<)8SKFU(I,*)*TE6P«*3 
?7*0  r-n(T,K)*SOCt,*l»TE"P»TE>*P 
2790  C.iK  riwuE 

l-n  3*0 0  1*1 . 12 
TF”P3AV(I.K) 

T‘FAsSO(I»k) 

T'l  J(3Dtl.'<)-TFMP*TePP/»6rN.'!)/(ANvflS*l,l 
IF (thp.lt.O.)  rnpeO. 

A  V ( T ,1) arEHP/ A6T»J 
Soil ,*)*THP»*.3 
T-PeSKE*!!.*) 

S*Ff  (l(«)"0. 

IF  (Hi  (1(0  .IE  .  ,nO05)  SO  TO  2800 


2207 

2208 
2209 
2210000 
2211 
2212 

2213 

2214 
2219 
2216 

2217 

2218 

2219 

2220 

2222 

2223 

2224 

2225 

2226 
2227 
2229 
2229 


2232 

2233 

2235 

2236 

2237 
2239 

2239 

2240 

2241 

2242 

2243 
224£ 
2203 

2246 

2247 
2249 
2249 

2230 
2250-1 

2231 

2232 
2253 

*  2254 
2253 
2256 
2237 
2236 

2259 

2260 
2261 
2262 

2263 

2264 

2265 

2266 
226T 
2268 

2269 

2270 

2271 

2272 

2273 

2274 

2275 

2276 
22T7 


EtHlPtT  6 
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:>k£W(I,K)»(*nyRS«*2»thP-3.«AFYRS*T£MP*TnPAt2,*TE.'1P»*3) 

1/f*NY9S*£iNY:(3«l,)«(»hY»5*?»)*  30{I,K)*»3) 
aaoc  cn»Tj»aie 

sate  continue 

LOTAP»N0TAP 

TTRNSal 

IF (IRC0N,LE.O)  GO  TO  2930 
IF(NCSTT.GT.O)  F,n  TO  1 2<»0 

C  Pfll'lT  *DJU»TED  FREQUENCY  STATISTICS 

282C  iPrTg(6.lO) 

'•‘RlTE(6.2*30) 

JOXC  FPPM*T(/30M  YOJI'STEO  FWECUENCY  STATISTIC*) 

M>tTF  (6, *90)  (NO(I),I«1 ,1?) 

on  PMO  K»».ST**,NSTA 

mpITE  (6,1010)  IjTACK), (Av(I,K),I«l,12) 

8  R I T  E  (6,1020)  (sn(l,K), loi.iaj 
irite  (6,  to3n)  (sxt^d, ‘<),lal,l2) 

■>o tte  (6,  loan)  (ond.K),  I«i  ,12) 

2s«c  CoftimiE 

c  PNImT  CUN5ISTEM  CCRRElATinN  matrix 

! TRUSS  1 
r.n  to  2170 

20Fo  IF  (IPCHS.UE.O)  GO  TO  2*70 
C  PtlNCH  FPEGUEACY  STATISTICS 

DP  2860  KSHSTAX.KST* 

‘•OITE(7,80)I6T*(K),  CAY  (1,8)  ,  t»l,l2) 

*  R  1 T  F (7,80) IsTA(K) , (SC (I, A), 1*1,12) 

hBlTF  (7,S0)  I:;T»(K) ,  (5XF.A  (I,K)  ,I3»,I2) 

RITE (7,90) 1ST  A (K), <00(1, K), 101,12) 

2A6 C  CONTINUE 

c  CO^PIOE  COMBINATION  FLONS 

c 

287c  IF(MCPMR.LC.O)  GO  TO  2910 
lTMPsl2«|*Yr,3»l 
0"  2900  M=2,ITmP 
0”  2090  K  XT  l  ,  ‘1C  OHO 
<SF  YtMSTA 

I  T  r  3  N  S  T  A  C  (KX,  IPASS) 

)(M.K)aO, 

Cn  2P«0  L°1,TTP 
irFPRexSTAClTV.L.IPASS) 

2  68C  Ifx.A  )aO(H,K)*r.(fl,ITtFP)*CSTAC(KX,U,IPAS*) 

209C  CONTINUE 
29CC  COI.TINue 
c 

CO  •  ♦  •  *  »  •  H*X  AM)  MIN  RECC6STITUTE0  FUCHS  •  »**«•* 

2910  oso 

T  T  f>  N  $  3  1 

JR  (MXRCS  ,LE  »0)  CC  TO  2*170 
IT, Rs\y9S 

2920  IF fTTHP.LE.O)  00  TO  2930 

-,T\*1 

Vls«x»C9 

lTVp5jTHp«MY«C$ 

IF(UHP.GE.O)  go  TO  3T3C 
Itmpsmx9C3»IT8F 
*■  t=TTMP 
I T“O*0 

GO  TM  3T3O 

293c  lF(IGMil.NF.?)GO  TO  302C 


C  C0MPOTF  OENEPAumO  STATISTICS* 

xW  T  T£ (6, 130) 

On  3000  K»1,nSTA 

C  AyFi;  AGE  CORPEU*  T  ICN  COEFFICIENT 


On  2950  L»1,K 
UT3L*K3TA 

8»vrx,u)*o, 

on  29*0  I»l,l2 

T"PSRA(J,K,L) 

IFfL.5E.X)TMp»RA(!,X,UX) 
29«C  °*v(K.L)«RAV(K,U)*TMP 
Kav(K,l)»RAV(<,L)/12, 


2278 

2279 
22*0 
22*1 
2282 
2283 
22*9  » 
22*3 
22*6 
22*7 
22*8 
22*9 

2290 

2291 

2292 

2293 
2299 

2295 

2296 

2297 

2298 

2299 

2300 

2301 

2302 

2303 
2309 

2305 

2306 

2307 

2308 

2309 

2310 

2311 

2312 

2313 
2319 

2316 

2317 


2320 

2321 

2322 

2323 
2329 

2325 

2326 

2327 

2328 

2329 

2330 

2331 

2332 

2333 
23J9 

2335 

2336 

2337 

2338 

2339 
2J«0 
23«t 

2392 

239J 

2399 

21*5 

2396 

23*T 

2398 

23*9 
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Exhibit  6 


woJTE{6.?0)rsT*{K)»ISTAa),R*V{K,L) 

2Vic  continue 

c  AVEUAGE  LOG.'!  FOR  MET  aNO  DRV  SEASONS 

AvPXHO«AV(ll*K)TAvCl2,X)*AVfl,K) 
t«y (K)»J 
AvNN(K)a4VMX(K) 

IMN fK) Bl 

TMP»Av(12,K)+*V(l,K)+*Vt2,*) 

If (AVHX(K).GF.TMp)GO  TO  2160 
AVMT(K)aTMR 
IP*  (Klo2 
on  TO  2970 
?06C  A  v *N  (K ) 8  TNP 
IMN  (K)»2 

c  Asp  4VEHA0E  STAfOARO  DEVIATION 

2rT0  S0AV(K)3sntl,K)*SD(2,K) 

on  29R.J  1*3*1? 

SI  AW (K)BSOAV { K ) >S0 (I,K) 

TMF*AV(t-?,K)+4V(I-l,V)+AV(I,K) 

If (AVMXOO  .GF.THPIGO  TO  2980 
AvNX  CK)«THP 
tMX (  K)  ■  I 

29*)C  If  (AvmnIn)  .U.TMP)GO  TO  2990 
AVMN(K)sTMp 
IMS IK)»I 
?RRC  CONTINUE 

Avm*  (K)=AVMT  (K)/J. 

AVS*N(K)SAVMNifK)/3. 
jOAV (K) 330*V(f) /!?, 

tore  msriNiiE 

'RtTEtJ..l«01 

Dm  3010  X°1**(ST4 

iTPajMxfiO 

ITVOolMU(K) 

lots  ~CITF.<6. 120)  lSTA(lO  ,AVHX(K),  AVMN(k)  *S0AV(K)  ,M0(lTP)  *M0(ITMP> 

C  »»*««»  At*Pi  T  GENERALIZED  STATISTICS*  ***»*•*••*•** 
10?C  If  CtGNKL.I.E.OlSn  TO  3100 
On  TOGO  K  =  1.N*1TA 
<ti-W*NSTA 

C  InTePeEPIATE  MCMHS 

;-»>!PBlMv(i<)-rHX(K)«3 

Jf  (N'M*rtN.LT.0)N*«XMM«.NMXMN*l2 

NNNH]fSfc*HM)(MN 

On  30«0  lBi *  i2 

c  STANDARD  DEVIATION  UNIFORM,  SKE*  ZERO 

AK?WtI,K)=0. 
no f I. K 1*0. 

■>n(I,N)*SOAV(KJ 
on  1030  L*1,iSTA 

C  7ERn  CORRELATION  MITH  OTHER  STATIONS  AND  PRECEDING  MONTH 

lyiL*NsTA 
' A  I  T , M , Ln 1 aD . 

If (L.GE.M)GD  TO  3030 

c  Uniform  SERIAL  CCPRfL  INTERMEDIATE  MONTHS  And  INTER*STA 

LA  (I,K,L)sR*v<i<,L) 
f A (T,l,K)SRA(t,K,L) 

303S  CONTINUE 

~:A  rr,*,KTl3rlAV(KfK) 

■'»  (  T  .  N  .  K  )  a  l  . 

tn°c  continue 

c  Mf  j.,  tNp  septal  CLRREl*  met  and  DRV  SEASONS 

TMPBRAvlK.K) n.IR 
r f»  p*t  mP«  ,  3 
If (TMP.GT..RR5T  Mps.96 
If (TEHP,LT.O)TEMP«0, 

ITP*IM*(«I 

’•V(;TP,K)*AV!iX(K)»,t 

(ITF.M.MxJaTEMP 
Itp«ImxIk)-1 
IF  CtTF.lT.lt  iTpaU 
AW(TTP,F)*AVMX{K)*,J 
'AUTF.M.MXIsTEMP 


2390 

2391 

2392 

2393 
2399 
2393 
2396 

2337 

2338 

2399 
2360 
236t 

2362 

2363 

2364 
2363 

2366 

2367 

2368 

2369 

2370 

2371 

2372 

2373 

2374 
2373 

2376 

2377 

2378 

2379 
23*0 
23*1 
23*2 

2383 

2384 

2385 

2386 

2387 
2348 

2389 

2390 

2391 

2392 

2393 

2394 

2395 

2396 

2397 

2398 
2394 

2400 

2401 

2402 

2403 

2404 

2405 

2406 

2407 

2408 

2409 
2»10 

2411 

2412 

2413 

2414 
2*13 

2416 

2417 

2418 

2419 

2420 
2*21 


EXHIBIT  6 


•  20* 


2622 

If (TTP.LT.l) ITP=ITP»12 

2623 

Av(TTP,K)a*V'l*(K)».l 

2fl2« 

^ACITP,K,kX)aT£"p 

2423 

ITPsIMN(K) 

2426 

*V{TTP,K)3AVSN{K) 

2427 

9A(ITP,K,KX)aTpP 

242  e 

lTPalMN(K)-l 

2429 

If (ITP.LT.l) ITP»12 

2430 

Av(!TP,K)*AVmN(K) 

2431 

<5A(TTP,K,KX)sTHP 

2432 

ITPalNN(K)-2 

2433 

If  (TTP.LT.l)  lTpalTp.»12 

2614 

AVCrTP.f  JbaVnN(K) 

2433 

f.4(TTP.X,KX)aTl*P 

2436 

c 

M£i‘iS  FU«  HCfTHS  FPLL0HlNG  *ET  SEASON 

2437 

IF (NMXMM.LT . 1 ) GO  TO  30bfl 

2436 

Itpsimx(k) 

2439 

rrnp=NHXPN* i 

2440 

TFHPa(AVHX(K)«.l-»VMNfK))/TEMP 

2441 

01  1050  IX31.NPXMN 

2442 

THP3IX 

2443 

IaIMx(K)*IX 

2444 

IPII.GT.12) 1=1-12 

2445 

H)3C 

AV(I,N)SAVtITP,K).TEHP*TPP 

2446 

C 

«FA'jS  Fp9  MQk  TPS  FOLLOWING  QNT  SEASON 

2447 

*060 

IF  (MHPHX.LT. 1 )  CO  TO  309C 

2448 

I  TPc INu  CK) 

2449 

TFNP=NMNNXn1 

2450 

TFHPetAVMX(Kl-.l-AVMN(K)l/TEHP 

2451 

Pfl  3070  IXal.NPNPX 

2432 

T  MP* I X 

2453 

i=r»N(K)*rx 

2454 

U  (I. GT. 12)1  =  1-12 

2455 

707C 

Av(T.K)SAP(lTP.K)TTEMP«TNP 

2456 

108C 

CONTINUE 

2437 

lone 

IPNPLSO 

2438 

I PCON*0 

2499 

00  TCI  1730 

2460 

Hoc 

lFCNyOR.LF.O.AN0.NPP0J.LE.O.AN0.NPA3S.LE,l)  GO  TO  150 

2461  « 

CP  • 

a  *  2462 

n;\up3nsta 

2463 

,VA0:iNSTA+1 

2464 

On  3200  1=1. 12 

2465 

IPSl-1 

2466 

tf  dP.LT.l)  IP  =  12 

2467 

on  3190  k*1.‘I,9TA 

2468 

On  3 1 a0  Lal»o3T* 

2469 

c 

CONftfLAT IONS  IN  CUBHENT  nONTH 

2470 

IF  fL.GE.<)  f.D  TO  3120 

2471 

.:fl.,NVAB)  a  03L£(RAd,K,L)) 

2472000 

np  1110  L*=L.nSTA 

2473 

Lv  =  I.  a*‘(3T* 

2474 

IF  (LA.LT.O  9  (L  ,  L  A  )  a  RBLE(hA(I,1,LA)) 

2475000 

if  (LA.GF.a)  9(1, LA)  a  ObLt  (F  * ( i , l , LX ) ) 

2476000 

31 1C 

«  (LA ,L)a9(L.L*J 

2477 

r.n  rn  3l«o 

2478 

C 

CP9 vCL ATI I)Ng  NI7N  P PECEOINS  nOnIH 

2479 

3  l  ?C 

L«*L*n.9TA 

2480 

HL.NVA9)  3  n«LE(BA(I,a,L*)) 

2441000 

nn  Mlo  lAsl.nsta 

2482 

910.1*)  =  P«LE(B»(Ip»l.,LA)) 

2483000 

11  JC 

<(L*.L)aB(L.L*) 

248  4 

11  «C 

CONTINUE 

2485 

C 

>111313311 

2486 

CALL  Cf)OliT(H) 

2487 

C 

cnisasiii 

2488 

On  lllo  l*1»nsT* 

2489 

me 

■ft* (I.a.L)-H(L) 

2490 

IF (OTBMC .LE . 1 .)  50  TO  3170 

249| 

VKTTf.  (6,3160)  I, K.OTRNC 

2492 

3U-C 

FrtBi’ *T  (3«H  INCONSISTENT  CC'BREl  BATBIX  FOB  Ia,l3,«H  la, 12, 

2493000 
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m<  orpHs«.F6.3) 
itsns»2 
.-;n  ro  i?30 

11?c  IFJOTSHC.fiE.O.)  r,0  TO  31*0 
M.ITE{6,70)!,k,Dtrhc 
0TW«C»0. 

71*0  ‘'LCPTa.K)3(1..0TflHC)««,3 
HI <1 C  CONTINUE 
32PC  CONTINUE 

C  «*«««•  GENERATE  FtO*S  «*«»♦*•*****••»*»» 

If  (NPASS.t.E.l)  GO  TO  J2«o 
32 1 C  If  (ISTAT.Eq.nST‘T)  GO  Tf  3220 
<f»n  ( 1 3 T A T ) 

LntATSLSTAT+1 
f..i  m  32jo 

72? C  “BITf;(IRTAT)!l7T3)l,NSTA,  (1RTAJK)  ,K«J,N3TA) 

•  .IT  A  T  a  M  S  T  A  T  ♦  1 
I.ST*T»N.1T»T 

On  1230  K=l.NST* 

<nITEfl3rATl  ISTA(K),(AV(I,!0,sn(I,Kl»3KtwII#K)»DG<I,K), 

1  (3ET*n.»<«t)»L»l»NSTA),ALCfT(I.K)*l»l,12) 

3230  voTAT»nSTAT*1 
l  .1 T  A  T®  NS  TA  T 

Tf  (TEAS3.LT  . |<P A 3 3 3  GO  TP  200 
52«C  Mai 

Ioass®! 

>.tQ 

"A  AO 

IF  (NPqnj.uE.O)  GO  TO  3310 

CO  •  A  .  •  A  A  Psn.jCc  Tgo  fucw  SEQUENCES  *a»a»»»aa«a«» 
32*C  JAs TrRpJ.IVNA 

‘  !=LYBPJ-tY9A 
l  r  m  P  a  0 

I  TPpmTHPJ«IMnTp«1 

If(ITP.NE.O)  oc  TO  32*0 
Itvps12 

72fcC  ! F  (irp.LT.I)  ITPSITP+12 

OAs  (JA-n  ♦  t2*I  TP  +  1  «IT*P 
!;o  7290  K  =  t  .m'TA 

If  fsr'ITTP.'O  .EO.O..C°.',A.eo.l)  GO  TO  32»0 
f  ►  -I's  A  tnc  0)  *(>4(I  f  P,k)  )  ».4  34,>945 

•P  •  Ev  CK.)  a  (  Tf  •'I’.lVflTP.K) )  /SOT  I  TP,  X) 

;r  fS*EW(ITP,Kl.r0,0.)  GO  TO  32PO 
TfPPa,5ASl<E«(ITP,K)A0PNfV(A)+l. 

T  .Jct  . 

tt  (TFHO.GE.O.)  GO  TO  3270 

rt -ps  t-tepp) 

T>'«a  («7Hp) 

327 r  rv {K)c6.»  t Tf P.TEMpA*  n ,/3.)«l .)/*KEn(ITP»x)*8KEn(ITP,K)/6, 
:-n  t o  32<)0 
3  ?  4  C  : p  f.  v  ( K  )  a  0  , 

32Gc  Continue 

jxsrvppj-l 

c  NS  GE'IuEnCE  NT.,  H  a  MONTH  NO,*  JX  ■  7E*S  NC, 

730C  'iaf*t 

.»0  TO  3330 

C  3  T  *  f T  a  I TM  2Epr  OfVIATlOs  AT  ALL  STATIONS 

331C  '  "  1320  K  *  t  » K  .7  T  A 
732C  '.'Pi-EvOOsO, 

c  generate  2  rfAOS  FOB  OIJCASOING 

’  I  -t? 

I Y  s  •  2 

333C  If  i>'PA39.LE.l)  C:  TO  3uCC 
If (1PASS.GT.13  r  J  TO  3J0C 
■'F..IM-  ISTAT 
OTAP»0 
t  <U  Apao 

33®C  Pf  •  I  no  IQTAP 
LOT ApaO 

''F.10(I3TAT)NSTY*,N3TA,(ISTA(kJ,K«1,N3TA) 

•  s  r  ysnoTx *»  1 

If (NSTX.lE.O)  GO  TO  3380 
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I  T»”16  J*  1  2*  1 
in  3.370  K=1»HATX 
If CIPOSS.Lf.t)  GO  TO  33*0 
338C  «P*0(TnTAP)  JTEmp,  (0(t<,K)  ,Fa?,lTP) 

L0TAP«li5T*P*1 

IP(ITPMP.NE.ISTA(K))  GU  TO  3330 

i3*c  t.r *o ( ist at)  jp*  (av(t»>o .sodiKj.SKewdxKNoajijX) ,  (bet*  t i » k « u >  #l»i 
I.MSTA).ALCFT(I,K),|'1»1«) 

137C  COaTINuE 
33*0  An  338p  K=»pATXX,NSTA 
IST«P«I9TAP*1 

tPCN.GT.O)  1pPE''fX)«03TAp{I9TAP) 

3  31C  AOCISTAT)  IP,  (AVCI,K),sn(I.K).3KEw(I»K).0Q(I,K),(HETA(I,X,L)#iai 

)  .-JSTA)  ,  A(.CFT{T,K)*IS»«  '-?) 

CB  *****  *  f,  E  N  F  0  *  T  F  cohsflateo  standaso  oetiate  ********** 
3«O0  IP(rPASJ.E'J.I)  JX  T  MPa  J  X 
.,rO*i h»hCCm"  (  jPaSI) 

‘.TEOhciTnO'iC  I P  a  3  3 ) 
in  3420  K»l.‘»9TA 
in  xato  1=1. 12 

tVG (I.K)aO. 

GOVfl.KlsO. 
xiilC  Continue 
3u?C  CONTINUE 

IffP.I.E^O)  GO  TO  loao 
ppITE(b.iO) 

.HITE  (6,34iO)  N 

3'j3C  rmiMAT  (2?h  <;FnEPATED  FLOWS  POP  PF3I0D.I3) 

TP(VPAS8.GT,1)  witITE(6,2630)  IPA33 
3442  on  1510  J  = 

tal2«tJ-l)»l 
on  ViOO  1*1,12 
•  =  -  +  1 

IE  (NSTx.tE.O)  Gfl  TO  3460 
Of,  3«3o  x  a  i  ,  9  T  X 
3„5r  ;psPv(K)afl(M,K) 

(Uf-C  IP  fH'. V.E  00  TO  35G0 

in  xapo  ksiigtxx »n3TA 
c  HANonn  COMPONENT 

TEEPaO . 
in  3370  L  = 1 »  * 

TrEParFwP*PMC.FN(IXX) 

3  a  T  c  T  P*1  psTp,. a  ■  .>  Mr.fe  f.  {  ;  x  X  ) 

T  E  EPa  T  EMp  *  AL  C  P  T CI,S) 

in  «a«0  L B 1 »  m 5 T A 

ld«n  te  >  (>ST6“P*‘  ‘  T1  ( I .  «  »l)  *IJP"E  V  (L) 

‘•Ki'ifl  .K)aA\;G(J»K)*TErtP 
■u '  v  ( I  .K)  a  Six  ( t .  K)  +  TEMr»Tt»*P 
"  f H  ,  x ) sTEHP 
'■  pip  V  (X  )  aTE'lP 
ThOC  r  •>.  TIMJE 
3snr  C-.  TTNUE 
Tr,  I  C  CONTINUE 

Ie  (NPA.99.lE.l3  GO  TO  3530 
35:>C  IP  (LOTAB.Eo.m'JTAP)  GO  TC  3530 
• f »n  (IOTAP) 
l.oTAPaLOTAp*! 
in  t n  3520 
353C  itps>.J*12»| 

[-.T  apbHT4P«*iSTA*NSTX 
in  1540  K  =  '.8TX  1 ,  ►  8T* 

'I  ft (IOTAP) IST*(K) , (0(P,K),na2,ITP) 

■■  IT  APBHQT  A P ♦  I 
I  -I  TAPajSTAp*! 

IE  n9T*P.r,T.K3TAP)  00  TC  160 
35«n  sT»P(l3T»p)=rj(lxP|X) 

355C  AM.nr.SNj.jA*l 

in  1fcT0  **> 9TXX.N8TA 

IE  (5J*JXTMP.r,r,0)  HQ  I  T  E  (  6 . 2660)  (MO  ( I ) »  !•* »  *«* 
rn  33*0  I « 1  ,  i 2 

Avr,fi,K)sAvr.(I,X)/ANLCO  ,rc-  « 

5nvfI.xT»((3,,v(l,i<)«AvG(I,P)«aJ»AFL0G)/»NV.C<»)a»»5 
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3560  fnsftNU* 

.Iva.fxT“;J 

on  *660  joja.nj 
JxaJX#t 
Fal2*J«l t 

IF  ( JX.IE.O )  00  TO  3660 
I  TP*0 

on  1630  1=1*12 


c 


(«M*  1 

TP  00  TO  3600  u.DT*»c 

TRANSFORM  TO  LCG  PF19306  TT Ft  111  VABJATE  (FLOW} 
T“PBSKE«(I,t<) 

tF(ANL0fi.OT.l9..AN0.3Dv<r.<'  J  .St. 0.1 
J  «fM*K)»fO(N.K)«AVG(I,K))/80''CI#tO 
IF  (XFp.EO.O.)  GC  TO  16CC 


T>«Pa(  (TNP«(ncN,K)»THP/6.  j/6,4t ,  J  «*3  •  i  « J  *2* /THP 
T£t*P3f.?.)/.9xEI>(I,jC) 
tF(SKF.W(T»K^  35A0,  3600,3590 
tF  (THP  .6t  •TEmP)  ?>'PbTEnp 

r-n  TO  3610 

3S<»C  IF  (  TMP  *1.  T  ■ TE “0 )  THPoTEhP 
on  to  3610 

?60C  T“P»a(H,X>  „  .  .  ......  w» 

16  (C  tF(T,xP.GT.3..AK0.S0a,K).GT,.3)  T»P»2  ,  ♦  1 1  MP1^  t )  *  .3/90 1 1  *  *  > 
TH?aTtM»*90{t,K)#AV(I*K) 

3.**T«P-oa(I.B) 


?6ie 
26»0 
26#  1 
?6#2 
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26#7 
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2631 

KlTWOBE*  2652 
265J 
265# 
2659 
2656 
265? 
265# 

2659 

2660 
266! 
2662 
2663 


JTFPCKStJf  (Kt  IPA33) 

IF (tTHF.LE.O)  CO  TO  3630 
TFXt>«C. 

an  1620  Lal.rT,,p 

Lv=ISTO<.>..  in*##) 

3h?C  tF“PSTEmptq(x,(,XJ 

CF(n(H,K) .lT.Temp)  n(xtKl«TExP 

7ft  JO  IFfn(M,K).LT.n,.ANP,QXIK(T,Kl.GE.0.1  QlM«K)*0, 
T6#C  t"f T)«OfX,K) *.5 
tTP«tTf>*IO(l) 
tiS«C  CONTINUE 
C 

in (1 31  at  TP 

--ITF  16,100)  1ST* C K 3 , jx, ( J07I) *I»1* 13) 

IF ( fPCHO.LF.OlGO  TO  3460 

,c fTff7 ,2730)  ISTA(K),  J X  ,  ( TQ 1 1 1 , 1» 1 / 1 2> 

1640  msTINUE 
367c  CONTINUE 

IF  (xernO.t.E.n)  r'U  TO  3720 
!>n  5710 

'« s  1 2  •  J  •  1 1 
nn  3700  1*1 .  12 

l'3X*l 

c  COMPUTE  C0NUT6»TIOn  FtOx.T 

'in  3660  KVal.MCO."# 

►.  s<  y  ♦  n  3  T  a 
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ITPsNJTAClKX*  IP*33) 

(«,ObO, 

■  n  TMo  i, s l ,  I  TP 
ITF“PaK4TAC(K A*L,TPASS) 

k)  Brt(  '»,‘<l-*Ofx,lTePF)«Cr.T»C(KX,U,  124.13) 

continue 

C^kTIKoF 
Onr,  rt*uF 

IF{N,tT.S"#nj!  50  TO  3250 
IF  (ntmxG.I  E.  >1C0  to  J#aO 

*  ,  .  «  •  M&X  #60  “IN  C66£8ATf  0  F(.0“8  **•#*»«* 
IF ( JX .t_F .0 )G  I  TO  1870 

SKIP  f*x“!N  IF  #E“*  TNJNG  TEA9S  INSUFFICIENT 
rF(JX.Gt.0.*,N0.NJ,LT.6T“XC)60  TP  150 


iTf.'.SBl) 

ITFbNSTAmMCOf.9 

on  .1800  x»n?txx,itp 

Fit  0#ce4 0*9  «C  1“12,  MAX  NO  13f 

on  1? to  r*i»t5 


6*N0  1#»  3*“N0  15 
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17«C  V<3fI#K)"«T 

C  HIM  CALENDAR  MC  14-27,  MIN  «0  2«,  h»MO  £9,  *a,MU  JO 

i'n  3730  1316,30  *  3  "u  JO 

375c  3mO(T,K)bT 

C  TMP  8  6. MO,  TEMP  a  94. MO  VOLUME,  THPA  »  |.Mn 

TEmobo.  ’ 

TMPaO. 

*VGO(K)30. 

soafi 

"31 

IF(ITRNS.GT.O)  Hs(H-1)«mxHC8«12*1 
')0  3730  J»1,MJ 
on  3780  I»lfl2 
IJ»I»13 

"3M*1 

T  M  P  A  a  fl (  M  ,  K  ) 

'VG0CK)8*VG0(K)»TMP» 

Mfl3M0*l 

IP(THPA.GT,8HfltI,K)13MQ(i,K)sTMP* 

IF  <TMPA.LT.9Mn(IX,K))3Mc(!x,K)sTMPA 

IP  tTMPA,r.T.3n0(t3,K))8M<J(13,K)sTMpA 

IP{TMPA,LT.3Mll(2P,K))8nO{2»,K)aTHPA 

tmp3tmp*tmpa 

TEmPsTENP+THMA 

IP(''.IT.8)GQ  70  3760 

T'lpaT-r-llfH-M,*) 

tP(rHP.Gi.3MT(in,K))SMn(iA,K)sTMP 
IP (TMP.UT.9H  Of 29, Ml >3*0  ( i?  9,  K  jay  HP 
IP  (U.LT.56)GC1  ro  1770 
rp-<psTEMp«a('<*5M,K) 

IP  ftEMp.GT,SM'J(13,K))SftC  (15,K)aTfcNP 
IP  (TpMp.LT.8-(3  (10,  K)  J3MH  (30,K)«TEMP 
•'.n  Til  17  AO 
3760  iM<Hl«,K)»TMP 
T77c  -juij  f  1  5  ,  K  1  BTEiP 
3 7 AC  OMstimuE 
J75c  continue 

c  average  monthly  flow 

TP  :'P3Sfl 

»VGO{K)8AVG0{l<)/TEHP 
3mOC  rNNTtV!'jE 

<“[TF  (4.10) 

IP  (ITPNS.Gf  .'))*9IT£(b,3?10)N,NJ 
3MC  r-PMAT  (/27H  “ATImijm  vuhj*F3  Fc«  PE«IOO,13,3H  OP, la, 

1  4PM  YEAK3  13F  MtCUMOEO  AKfi  PEC  ON  8  T I T  UT  EO  PLO«S) 
IP(TrPM3.LE.O)MHITE<6,3P^O)N,NJ 
1MC  ^n^MAT  (/2/M  *A»IMUM  VOLLMtS  FcR  PER  103, 1 3 , 3h  OP,I«, 
l^M  YFAHS  OF  SYNTHETIC  FLFm3) 

,0  I  TE  (6, FIO)  (‘10(1).  1*1,12) 

I  fFasSTI  +  MCO’’" 
in  1440  K8M3TYT, ITP 
I tf  xpaAVGO (K) ♦ ,5 
1 n  inJO  1*1, |0 
103C  TnCT)3SM,i(I#K)Y,5 

»"t  TE (6,840)  I M  A  (K),  (lQfI),T«l,15),  ITE«P 
J84C  C " *  T INuE 

TJ  fK (4,650) 

•PITE (6,610) ("0 (I), 1*1, 12) 

"0  1440  K6N.6TYV.ITP 
00  <«5o  131,15 
163C  inf  l)68M0(I*l5,K)a,5 

nPITP(6,840)l3TA(K) , (IU(I),l9i,13) 

38*C  CC'-TIMuE 

C  TRANSFER  t'ACK  TO  REC0N3TITUTF.0  210*8 

IF<TTRn3.GT.0)GQ  TO  2920 
1«7r  NT  a  NYHXG 
SO  TO  3990 
3680  ‘J  a  KyR 

3o°n  ;f(mpajs.le.i)  go  to  J9oo 

'P»*«*IPA81Tl 

IF(M.EO.O,AN').IPA88,Ll,APASS)  g(J  JQ  jug 
IF (TPAjs.LC.mpass)  go  TC  334q 
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IPA93»t  2781 

C  Gu  rO  Mth  JO1’  27*2 

3<7<V  IE(NYRli.l.E.O)  00  TO  150  27,3 

IE(RJ.GT.MrRf,)NJ=NYHG  2™ 

VPS3NYRG-»:  J  2785 

on  TO  3100  2788 

End  2787 

sMBonume  chouhrx)  ' oo i 

OTKENSlOs  Htli>),<»Cin,in,RX(10,n)  10'12 

DOUBLE  PHFCI3ICN  R,8,Rx  1093 

COMMON  0TR.iC,NINDP*8  1004 

*.«ARONINOP*1  1005 

On  20  J3 1 1  s*  I n 0 p  1006 

1007 

On  to  K  =  t.NWAP  1008 

tC  nfJ.OsOxfj.K)  1009 

SC  COM INUE  1010 

KfRINDP.GT.DSO  T9  30  1011 

'(11=9(1,71^1(1,1)  1012 

Ofi’MC««(l)  *0(1)  1013 

RETURN  101# 

c  *».«**  nfc-«iVEl>  MATRIX  A********************  1015 
30  i)0  #0  K'2.'JV*M  1016 

#n  J fl ,K)««(t ,*) /K(l , l)  1017 

,7,1  no  K  =  2  • ‘it  |0P  1018 

t TB=K»1  101# 

nn  60  J*K,‘.TuDP  1020 

r,ii  so  x=i»irD  1021 

U  rK  « I  1022 

5C  ,(J»*)*R(J,K)*'1(J.L)*R(L,K)  1023 

JE  (J.EQ.K)  Su  TO  60  102# 

;(K,J)aP(.l.«)/R(K»K)  1025 

*C.  CnNTtMlE  1026 

nn  70  1=1, tTP  *  1027 

L-nX-I  1028 

7n  -  f«  VAH=  I  (k,NVAB)-R(U,6»*P)*R(K,L)  1029 

f5«P*0A8Sfit(H«Kj)  1030 

jr  (TliHp.r, r ,  .oooor.l )  GP  T r  80  1031 

1T#»CB1.1  1032 

->:TliK6  1013 

cr.  i,;<,NVAP)3!?(k.NVA5)/M(K,K)  103# 
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EXHIBIT  7 


INHJT  DATA  723-X6-L2340 


CARD  VARIABLE 


COMMENTS 


A 


Three  title  cards,  first  must  have  A  in  column  1. 


B 


First  specification  card. 


C 


1.  IYRA 

2.  IMNTH 

3.  I ARAL 

k.  MXRCS 

5.  RYRG 

6.  NYMXG 

7.  NPASS 

8.  IPCOi 

9.  IPCHS 
10.  ROTA 

l.  1ICOMB 

2.  NTNDM 

3.  NCSTY 


-  Earliest  year  of  record  at  any  station. 

-  Calendar  month  number  of  first  month  of  water  year. 

-  Indicator,  positive  value  calls  for  statistical 
analysis  routines. 

-  Number  of  years  in  each  period  of  recorded  and  re¬ 
constituted  flows  for  vhlch  maximum  and  minimum  values 
are  to  be  obtained,  dimensioned  for  100. 

-  Total  number  of  years  of  hypothetical  flows  to  be 
generated . 

-  Number  of  years  in  each  period  of  generated  flows 
which  maximum  and  minimum  values  are  to  be  obtained, 
dimensioned  for  100. 

-  Number  of  consecutive  passes,  each  pass  consisting  of 
a  new  group  of  stations  vhlch  can  be  correlated  with 
specified  stations  in  previous  passes,  dimensioned 
for  5. 

-  Indicator,  positive  value  calls  for  writing  recorded 
and  reconstituted  flows  and  generated  flows  on  Tape  7. 

-  Indicator,  positive  value  calls  for  writing  stat¬ 
istics  on  Tape  7. 

-  Number  of  stations  at  which  flows  are  to  be  generated, 
not  required  if  flow  data  are  supplied.  NSTA  +  NCOMB 
(C-l)  dimensioned  for  10. 

Second  specification  card. 

-  Number  of  combinations  of  stations,  the  totals  of 
which  are  used  to  obtain  maximum  and  minimum  flows, 
dimensioned  far  2.  If  positive,  provide  D  and  E  cards. 

-  Number  of  tandem  situations,  compares  sum  of  monthly 
values  of  upstream  stations  with  downstream  station 

and  adjusts  if  value  Is  less  than  sun  and  that  station's 
value  has  been  estimated  or  generated,  dimension  for  10. 
If  positive,  provide  F  card. 

Number  of  consistency  tests.  Adjusts  standard  deviation 
of  a  dependent  station  in  tandem  with  an  independent 
station  to  prevent  frequency  curves  from  crossing, 
dimensioned  for  10.  If  positive,  provide  G  card. 
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CARD 


VARIABLE 
C  (Cont'd) 


4.  IORRL 


5.  HPROJ 

6.  IYRPJ 

7.  MTHPJ 

8.  LYRPJ 


D 

1.  If  STAC 

2.  I STAC 
E 

1.  NSTAC 

2.  COTAC 


F 

1.  I STB 

2.  IfSMX 

3.  ISTT 

G 

1.  ISTX 

2.  I STY 


H 


comments 


Indicator,  +  1  calls  for  reading  generalised  stat¬ 
istics  and  using  for  generation*  +  2  calls  for  com¬ 
puting  generalized  statistics  froa  flew  data  and  using 
for  generation. 

Humber  of  projections  of  future  flows  froa  present 
conditions*  usually  0. 

Year  of  start  of  each  projection. 

Calendar  month  of  start  of  each  projection. 

Last  year  of  each  projection,  number  of  recorded  and 
reconstituted  years  plus  number  of  projected  years 
dimensioned  for  100. 

Identification  of  combination,  IfCOMB  (C-l)  sets  of  D 
and  E  cards. 

Nuafcer  of  stations  in  this  combination,  dimensioned 
for  10. 

Station  number  (RSEAC  values). 

Combining  coefficients,  IfCOMB  (C-l)  sets  of  D  and 
E  cards. 


Same  as  D-l . 

Coefficient  of  flow  used  for  adding,  corresponds  to 
respective  items  in  D-2. 

Identification  of  tandem  situation,  NTNDM  (C-2)  cards. 

Station  number  of  downstream  station. 

Humber  of  upstream  stations,  dimensioned  for  10. 
Station  number  of  upstream  station  (WSMX  values). 

Identification  of  consistency  test,  HCSTY  (C-3)  cards. 

Independent  station  number. 

Dependent  station  number . 

Flow  data,  cards  in  any  order,  oalt  if  I  ANAL  (B-3)  is 
not  positive,  follow  all  flew  data  cards  by  1  blank 
card  (I  card). 


1.  Cols  2-4,  Station  number 

2.  Cols  5-8,  Year  number. 
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CARD 


VARIABLE 


COMMENTS 


H  (Cont'd) 


I 

J 


K 


3.  Cols  9-14,  15-20,  etc..  Flow  In  desired  units.  Units  should  be 
selected  so  generated  flows  will  not  exceed  999,999.  Use  -1  for 
missing  record.  If  record  for  entire  year  is  missing,  omit  card 
for  that  year. 


Card  blank  after  Col  1  to  indicate  end  of  flow  data, 
omit  if  IANAL  (B-3)  is  not  positive. 


1 .  NCOMB 

2 .  NTNDM 

3 .  NCSTY 

A.  NE7TX 

5.  ISTA 


Identification  of  stations  in  previous  passes  to  be 
used  in  current  pass,  supply  only  if  NPASS  (B-7)  is 
greater  than  1.  The  variables  NCOMB,  NTNDM,  and 
NCSTY  apply  to  the  current  pass  only. 

-  Number  of  combinations  of  stations,  the  totals  of 
which  are  used  to  obtain  maximum  and  minimum  flows, 
dimensioned  for  2.  If  positive,  provide  D  and  E  cards. 

-  Number  of  tandem  situations,  compares  sum  of  monthly 
values  of  upstream  stations  with  downstream  station 

and  adjusts  if  value  is  less  than  sum  and  that  station's 
value  has  been  estimated  or  generated,  dimension  for  10. 
If  positive,  provide  F  card. 

.  Number  of  consistency  tests.  Adjusts  standard  deviation 
of  a  dependent  station  in  tandem  with  an  independent 
station  to  prevent  frequency  curves  from  crossing, 
dimensioned  for  10.  If  positive,  provide  G  card. 

Number  of  stations  from  previous  passes  which  are  to 
be  used  with  the  additional  data  in  current  pass  as  a 
means  of  maintaining  consistent  flows  between  groups  of 
stations,  number  of  stations  from  previous  passes 
plus  number  of  new  stations  dimensioned  for  10. 

Station  number  of  station  In  a  previous  pass  which  is 
to  be  used  in  current  pass  (NSTX  values).  Must  be 
in  same  order  as  stations  first  appear. 


Note:  Flow  data  for  current  pass  supplied  ns  described  for  H  card  and 
folLow  data  with  a  blank  card  (i  card),  supply  NPA3S-1  sets 
of  J,  H,  and  I  cards  (also  D,E,F,  and  G,  if  necessary)  when 
NPASS  greater  than  1. 

Preceding-month  correlation  coefficients  for  first 
station,  omit  if  IANAL  (B-3)  is  positive  (NSIA  cards). 


1.  ISTA(K) 

2.  ISTA(L) 


-  Cols  2-U,  Number  of  first  station. 

-  Cols  5-8,  Number  of  station  from  1  to  NSTA  (B-10)  on 
successive  cards.  If  IGNRL  (C-**)  =  l,  only  first 
card  is  used. 


3.  RA(l,K,LX)  -  Cols  9-1^,  15-20,  etc.,  Correlation  coefficients  for 
successive  months  between  flows  at  first  station  and 
preceding-month  flows  at  stations  from  1  to  NSTA  (B-10) 
on  separate  cards.  If  IGNRL  (C-U)  =  1,  only  general¬ 
ized  coefficient  (in  cols  9-1*0  is  given. 
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CARD  VARIABLE  COMMENTS 


L*  Current-month  correlation  coefficients,  omit  if  IANAL 

(B-3)  is  positive,  (NSTA-l)  pairs  of  L  and  M  cards. 

1.  ISTA(K)  -  Cols  2-4,  Number  of  station,  progressing  from  K  =  2 

through  NSTA  (B-10)  stations  on  different  sets  of  L 
and  M  cards. 

2.  ISTA(l)  -  Cols  5-8,  Number  of  station,  progressing  on  different 

cards  through  all  stations  from  L  =  1  to  K-l. 

3.  RA(l,K,L)  -  Cols  9-14,  15-20,  etc..  Correlation  coefficient  for 

each  successive  calendar  month  between  flows  at  station 
K  and  concurrent  flows  at  station  L  (12  items).  If 
IGNRL  (C-4)  =  1,  only  generalized  coefficient  in  cols 
9-14  is  given. 

M*  Preceding-month  correlation  coefficients  for  remaining 

stations,  omit  if  IANAL  (B-3)  is  positive.  Faired  with 
L  card. 

1.  ISTA(K)  -  Cols  2-4,  Same  station  number  as  on  corresponding  L 

card  (L-l ) . 

2.  ISTA(l)  -  Cols  5-8,  Number  of  station,  progressing  in  same  order 

on  different  cards  through  all  stations  from  L  =  1  to 
NSTA  (B-10).  If  IGNRL  (C-4)  =*  1,  only  card  with  L  =  K 
is  used. 

3.  RA(l,K,LX)-  Cols  9-14,  15-20,  etc..  Correlation  coefficient  for  each 

successive  calendar  month  between  flows  at  station  K 
and  flows  in  preceding  month  at  station  L  (12  items). 

If  IGNRL  (C-4 )  =  1,  only  generalized  coefficient  in 
Cols  9-14  is  given. 

N  Generalized  frequency  statistics,  omit  if  IANAL  (B-3)  is 

positive  or  IGNRL  (C-4)  does  not  equal  1. 

1.  ISTA(K)  -  Cols  2-8,  Station  number  for  NSTA  (B-10)  stations  on 

successive  cards  in  same  order  as  supplied  by  L  cards 
(L-l). 

2.  AVMX(K)  -  Cols  9-14,  Average  mean  logarithm  for  wet  season  (3 

months ) . 

3.  AVMN(K)  -  Cols  15-20,  Average  mean  logarithm  for  dry  season  (3 

months ) . 

4.  SDAV(k)  -  Cols  21-26,  Average  standard  deviation  for  the  12  months. 


*  Sets  of  L  and  M  cards  are  required  for  each  station  from  K  ■  2  to  NSTA. 
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CARD 


VARIABLE 


COMMENTS 


N  (Cont  M ) 

5.  MOMX(k)  -  Calendar  number  of  last  month  of  wet  season. 

6.  MDMN(K)  -  Calendar  number  of  last  month  of  dry  season. 

0  Mean  logarithms,  omit  If  IANAL  (B-3)  is  positive  or  IGNRL 

(C-4 )  equals  1. 

1.  ISTA(K)  -  Same  as  (M-l). 

2.  AV(l,K)  -  Cols  9-14,  15-20,  etc..  Mean  logarithms  for  successive 

calendar  months . 

P  Standard  deviations,  omit  if  IANAL  (B-3)  is  positive  or 

IGNRL  (C-4)  equals  1. 

1.  ISTA(K)  -  Same  as  (M-l). 

2.  SD(l,K)  -  Cols  9-14,  15-20,  etc..  Standard  deviations  for  suc¬ 

cessive  calendar  months. 

Q  Skew  coefficients,  omit  if  IANAL  (B-3)  Is  positive  or 

IGNRL  (C-4)  equals  1. 

1.  IOTA(K)  -  Same  as  (M-l). 

2.  SKEW(l,K)  -  Cols  9-14,  15-20,  etc..  Skew  coefficients  for  suc¬ 

cessive  calendar  months. 

R  Flow  increments,  omit  if  IANAL  (B-3)  is  positive  or 

IGNRL  (C-4)  equals  1. 

1.  ISTA(K)  -  Same  as  (M-l). 

2,  DQ(l,K)  -  Cols  9-14,  15-20,  etc..  Flow  increments  for  successive 

calendar  months. 

Five  blank  cards  with  A  In  Col  1  of  first  should  follow  last  Job. 


Note:  Cards  K  through  R  are  not  required  if  cards  H  and  I  are  supplied.  Cards 
K  through  R  are  as  punched  by  computer  when  IPCHS  is  positive. 
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EXHIBIT  8  ' 

SUMMARY  OF  REQUIRED  CARDS 
723-X6-L23**0 


(l)  Supply  only  if  IAHAL  (B3)  Is  positive.  Repeat  u  card  for  each 
station-year  of  data  before  supplying  i  card. 
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SUMMARY  OF  REQUIRED  CARDS 
Continued 
723-X6-L23i*0 


L  designates  correlation  with  current  month  and  LX  with  preceding  month.  If 
IGNRL(C4)  =  1,  only  one  (generalized)  coefficient  is  given  following  station 
numbers  on  each  card  and  only  1  K  and  M  card  is  used  for  each  K  station,  with 
L  =  K.  Use  same  format  as  H  card . 

Repeat  set  of  L  and  M  cards  for  each  K  station  except  first. 

Omit  if  IAKAL  (B3 )  Is  positive. 
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